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We describe in this thesis an attempt to calculate ground 
state properties of dilute simple metallic alloys. The 
properties of interest are the ground state energy or the 
heat of solution, the equilibrium volume or density, the 
electron density distribution and charge ti*ansfer. These 
properties are of considerable fundamental and technological 
interest. Most of the recent theoretical work in alloys is 
concerned with the nature of electronic states (localized or 
itinerant) and their organization in energy (i,e, density of 
states). The methods developed (e.g. CPA, cluster expansion) 
are not suitable for calculating the properties mentioned 
first above. We apply here the general theory of the inho- 
mogeneous electron gas developed by Kohn and co— worksrs. 

In this theory, the ground state energy is shown to be a 
nn 1 qiifl functional of the electron densiiy , being a minimum 
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for the correct density. The theory has earlier been applied 
successfully to calculate work function and surface energy 
of metallic surfaces where the electron density falls to zsto 
typically in an atomic distance. In case of the alloy, the 
electron density changes involved are smaller. We calculate 
the ground state properties of an alloy in two ways. In the 
first, the electron density is parametrized suitably, and the 
ground state parameters are determined variationally, by 
minimizing the ground state energy. In the second, a 
Schrodinger like eq.uation derived from the variational prin~ 
oiple is solved to obtain wave functions and the 

electron density n(?). The basic physical idea behind both 
these methods is that of neutral, weakly coupled pseudo-atoms/ 
a metal atom A dissolved in B Carries with it, to a good 
approximation, its own electronic environment. Some density 
adjustment occurs at the A-B interface. 

In the first chapter of the thesis, we place in perspe- 
ctive the problem of calculating the heat of solution of 
alloys by summarizing the earlier largely (Qualitative ideas 
of size effect, electronegativity, etc, due to Hume-Hothoiy,_ 
Pauling and others. The application of pseudopotential 
perturbation theoiy is also reviewed. We then describe the 
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intLomogeiieous electron gas (density functional) theory of 
Kohn, Hohenherg and Sham, and outline the two major ways 
in which it can be applied to calculate the ground state 
energy of an inhomogeneous system consisting of an atom 
or ion A replacing the ion B in a lattice of B ions (the 
alloy is ^ with c<l< l) (Chapter 2). 

In the third chapter, we detail the variational method 
of finding the ground state energy (g.s.e,). The g.s.e# is 
assumed to be a functional of the electron density. All the 
terms contributing to the g.s.e,, namely kinetic energy, 
exchange, electron correlation, electrostatic interaction 
and eleotixjn ion interaction are included. The approxima- 
tions made are the following, 1) A local functional form 
is assumed for the first three terms. The forms and 
parameters are taken from the well developed theory of the 
interacting electron gas, 2) A term proportional to the 
sq.uare of the gradient of electron density is included. The 
coefficient of this term has been obtained recently in 
fairly sophisticated approximation schemes by various 
authors and by us, 3) A local pseudo potential approximates 
the effect of the ion on conduction electrons. We show how, 
given a density n(r), the energy can be expressed in a rapid' 
ly convergent form. The most important term involves an 



electrioally neutral neighbourhood of the ion h. Because of 
charge neutrality, other terms, -which can be calculated, are 
Small, Ve discuss how the heat of solution, lattice constant, 
bulk modulus etc, can be calcirlated. In a local theory with 
gradient terms the only two centre dependent term in energy 
arises f2X)m the electrostatic interaction. Since the inter- 
atomic terms contribute only about 10 percent of the local, 
density or volume dependent terms in simple metals, our 
approach of ooncentratiug on the latter is clearly reasona- 
ble, A more satisfactory and complete formalism is developed ' 
in Chapter 5« 

In this chapter (Chapter 3) the heat of solution has 
been calculated with two specific variational choices of 
n(r). Firstly, we vrork with one step model where the step 
location and height are variational parameters. One of these 
is fixed by charge neutrality, and we obtain analytical 
expressions for the heat of solution as a function of the 
single free parameter. The model is unrealistic in that 
because of the density discontinuity, the gradient term 
(approximation (2) of paragraph above) gives an infinite 
oontribution. However, comparisons -with a more realistic 
model shows that the step model gives the correct trend 
and order of magnitude of 41, In a second choice, we use a 



smooth, (tangent hyperbolic) distribution for nC^). This 
leads to a two free parameter model from which heat of 
Solution can be numerically obtained by minimization. 

The systematics of the heats of solution (partial 
molal enthalpy at infinite dilution) is discussed in Chapter 
4, for various atomic volumes and valence ratios of A and B 
(solute and solvent), .Ignoring the small interatomic term, 
we find that the heats of solution show well defined depen- 
dences on size difference and on the electronegativity diff- 
erence(which is defined here precisely in a new way). The 
numbers obtained are of the right size. We also calculate 
the charge transfer and attempt to make contact with various 
empirical and semi— smpiri cal descriptions of alloying para- 
meters, In a number of situations, the two parameter model 
does not work veiy well. The underlying physical causes are 
identified. 

In Chapter 5, we discuss the Kohn Sham Sch 2 x> dinger 
equation scheme for calculating nC^) and the g.s.e, Por 
pure metals, the scheme hsts been developed in a neutral 
pseudo-atom approximation model by Dagens, Ball and others. 

We present here a version of the neutral pseudo-atom model 
which is free from any arbitrariness, which clearly indicates 
the approximations involved^ and allows the corrections to 



be evaluated. The uethod is then applied to an 
alloy, and the g.s.e. can be interpreted as the sura of 
the energy of pseudo-atom A in the 3 matrix, the energies 
of the B pseudo-atoms, and' an interaction term. Preliminary 
results of gaxjund state energy for puro hi, pure Mg and H-Mg 
alloy are presented. 



CJhapter 1 


Introduction 

The role of electronic contribution to the free energy of 

alloys was recognised from the early ezperiments of Hume- 

fiothery and co-workers^. His observation led to an anpirical 

understanding of alloys in teims of size-difference# Yalency 

and electron per atom ratio. These famous rules based on size- 

factor, relative valence and the observation of specific 

crystal structures at definite electron concentration per atom 

are reviewed in the first section of this chapter in the 

light of their subsequent developments. This section also 

contains a discussion on the electro -chemical effect in 

2 

alloys starting from the Pauling's rule for estimating the 

» 

heat of mixing from the difference in eleotronegatfvities of 
the components. This rule although originally proposed for 
the ionic compounds, found a successful application in alloys^ 
in account of the recent attempt by Hodges and Stott to put 
this rule for the alloys on a theoretioal perspective has beai 
presented. This section, in essence, will tie in those as- 
pects of alloys which haVedeveloped mainly from the experimen- 
tal studies. 

The attempt to develop a theory of alloys from the first 
principles started from an urge to explain the striking 
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observations of the ear^.y experiments of Hur-.e-Eothery and 
others leading to Jones's"^ rigid hand model desoilhed briefly 
in Section 2, The shape of the observed f ermi-snrf ace lead 
one to question the correctness of the rigid tend model and 
refinements have been proposed since* In the fifties, 

5 

Friedel discussed the problaa of dilute alloys i,e, a 
single impurity in a metal, in terms of its electionic stru- 
cture, polarization due to screening, scattering phase shifts 
and virtual bound states (See also Mott^), ^ account of 
the essential conclusions of these investigations is presen- 
ted in Section 3« 

Prom here the development of the theory of alloys has 
taken two distinctly different paths. The first one starts 
from a simplified hypothetical model of disorder and studies 
the band structure, the density of states etc, of alloys with 
the help of approximations known as coherent potential appro- 
ximation (OPA) or the average t-matrix approximation (ATA) etc. 
These treatments although quite rigorous, have not yet been 
able to estimate such useful properties as the heat of 
solution, charge transfer etc, ¥e have completely left out 
this approach in this review. In the second approach tlie 
ground state energy (g.s.et), known to be a unique functio- 
nal of the electixjn density is directly evaluated. The 
evaluation uses the minimal property of the g,s,e, for the 
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cxjrrect de-.sity. Details of lj.e electronlo structure are not 
obtained in this scheme but it is specially suited for study- 
ing the g.s.e. (and hence the heat of solution), charge density 
and hence the charge transfer etc. It has been used successfully 
for metallic surfaces by Lang and Ebhn^ and others. In the 
limiting case of small periodic variation of charge density 
this approach coincides with the perturbation theoretic treat- 

7 

ment of alloys formulated by Inglesfield and others which 
has been presented in Section 4. 

1,1 Factors Affecting the Alloys ; 

The factors controlling the properties, of a given alloy 
like the heat of solution, crystal structure etc. can be broadly 
classified into two sets ; (a) the environmental factor, tem- 
perature and pressure which are associated with the entropy 
term and (b) geometrical and e_ectro chemical factors whose 
origin aro the electronic states in the alloys. The factors in 
the first category havo been studied thoroughly within the 
framework of thermodynamics and we exclude them from our 
study by considering the alloy system at T = 0°K, 

The second set of factors, electronic in origin, are at 
issue in this study and so in this section we will take them up 
for a review. The. last factor mentioned in this section namely 
the density mismatch effect, is not physical in origin but 
is an artifact of the theoretical construction of alloy. But 
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since a discussion on it is relevant to this study we have 
presented it as a separate factor. 

The factors discussed here are producteof experimental 
studies ahd so have been presented in different ways by different 
investigators. We will discuss them broadly under two cate- 
gories - the size factor and the electrochemical factor, 

1.1.1 Size factor effect : 

Hume-Eothery^ was the first to realise the effect of 
size difference on alloying which he expressed in his famous 
rule 'where the atomic diameter of the solvent and solute 

M 

differ by more than about 14 to 15 percent of the solvent, the 
size factor is unfavourable and the solid solution is veiy 
restricted. Atomic diameter is given by closest distaace of 
approach of the atoms in the crystal structure of elements'. 

This rule was arrived at while determining the limits of 
solid solubility of the elem«its of the B-subgroup with 
copper and silver. It was also observed that the lattice dis- 
tortion produced by an increasing size factor favours a less 
regular structure and so theliquidus and the solidus are 
depressed by the size factor effect. 

This rule, although negative in character, has been immen- 
sely successful in guiding metallurgists to search for systems 
with limited solubility where precipitation hardening can be 
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taken advantage for strength. jning. In the last fifty years 
this rule has not undeigone any significant change. Subse- 
(luent investigators have mainly concentrated on the definition 
of atomic diameter and a more reliable value for the limit 
in an empirical way. 

The theoretical estimate of the contribution of size 
difference to the free energy of alloys use the approxi- 
mation of an elastic continuum to represent the solvent and 
then considers the mechanical strain energy when the solute 
is inserted in this continuum. The model used is one in which 
a sphere of specified radius is forced into a hole of smaller or 
larger radius, the surfaces welded together, and the body 
allowed to relax in its self -stressed state. Bach solute atom 
is thus regarded as a centre of elastic deformation. 

This approach has met with a limited saccess in predi- 
cting solubility limits in several systems but has failed in 
most of the cases. The reason is not difficult to understand 
in view of the continuum assumption which breaks down on an 
atomic scale. Its success can also easily be understood if we 
expand the change of energy due to size mismatch, in terms of 
a Taylor series involving powers of the difference in size. 

It clearly tells us that if the size mismatch is small this 
approach may work. 
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Anoth r conceptual contr idiction we encounter with the 
continunm model is the long range nature of strain. But if one 
considers from the first principles as it will be evident from 
our study that all the terms involved in this estimate are 
only short range. 

Recently Hodges and Stott"^ have put forward a qualitative 
picture of size factor effect as distinct from electrochemical 
effect which the experimentalists have now^ to separate. In this 
study we have developed the idea quantitatively within the limi- 
ted framework of step model to be discussed in Chapter 2, 

Miedema and, co-workers in their empirical estimate of heat 
of solution in alloys have accountedfor size factor contribu- 
tion to energy by a teim which we fail to understand. In 
their first paper they assume that this contribution is pro- 
portional to the square of difference in electron density ( 6 n) 
of the constituents but in their subsequent papers they use 
a term proportional to (^n) 2/3 without any further justifica- 
tion. The method of computation of An has also defeated us 

because it finally leads to An = [(~r~a)^ - a) *]• 

■ ^^at /A - ^ at B' - 

In an attempt to understand the size factor contribution, 

we consider the neutral atoms as separately compressible 

entities in alloy. This can also be justified in view of 

Dagen’s^ finding that interaction eneigy between neutral 
* 

pseudoatoms is less than four percent of the energy of the 
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pseudoai/Oia for allcali metals. So for the dilute alloy 
case where the cell size of alloy is eq.ual to that of the 
solvent B we can write 


(Sa- 


3^E(r) 

ar^ Y = 


( 1 ) 


where and Rg are solute and solvent atomic radii respecti- 
vely and E is the energy of solute A with atomic radius r 
9 ^ B 

aiiti — ^ — is related to the bulk modulus of the solute A. The 

3r* 

validity of this simple ezpression has been examined in the more 
rigorous models under study, 

1,1,2 The electrochemical effect : 

fn contrast to the previous factor this effect has been 
subject to more widespread investigation, Hume-Rothery^ in 
his famous rules only observed the obvious fact that with the 
increase oi the electronegativity difference there will be 
an increasing tendency towards formation of compounds thereby 
restricting the solid solubility. The concept of electronegati- 
vity i,e, the power of an atom in a molecule to attract ele- 
ctrons to itself has been a guide to understand chonical rea- 
ctions but the precise meaning of this power in teims of first 

principles is still missing. In an attempt to give a consi- 

2 

stent mpirical scale of electronegativity Pauling in 1932 
correlated the difference in heat of formation of ionic com- 
poTinds with the square of the difference in electronegativity 
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of the constituent atoms. This when applied to the heat of 
Solution of alloys, was found to "be quite successful indica- 
ting the importance of electro-chemical effect in alloys. 
Since then several attempts have "been made to iinderstand the 
concept of electronegativity and identify it with some 
directly measurable property but each of them has met only 

Q 

With limited success, Miedema and co-worhers have identi- 
fied it with the work function and have tried to estimate the 

3 

electrochemical effect in alloys. But Hodges and Stott"^ on 
the basis of density functional formalism and some crude 
approximations have found the leading term in the express 
expression for the heat of solution to be proportional to the 
square of difference in the internal chemical potential of 
the constituents and hence identified the electronegativity 
with the chemical potential. 

If n (r) is the ground state electron distribution in 

s ^ 

alloys then the chemical potential of the systea is given by 

g 

\x as follows , 

I,', 

\x = — (2) 

3ng(^) 

' 4 - 

where E is the ground state eneigy of the system, 

s 

In equilibrium it is the same th 2 X)ugh out the system and so 
should be independent of r, Hodges and Stott assume that in 
the alloy A-B, atoms A and B are 'prepared’ (by appropriate 
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compressic i or dilation (from their equHibrium size in the 
pure state) so as to ’fit’ the alloy unit cell. Suppose the 
charge distribution now is rLj.(r). The electronic distribution 
will relax to the gixDund state distribution n (r). Prom fun- 

O 

ctional Taylor series expansion it can be shown that, upto 
second order 


B[n^] - Sg iff 

*“ * 


3^B[n] 

} 

3 n(r) 3nCi*) 


^ [n^CS-) - ng($)] 


[n^(r') - n^Cr*)] 
Taking derivative with respect to 




/ d^ { 


S^Bfnl ^ 
SnCi*’) 3n(r) 


- 0 

Substituting eqn. (3) iii (2) 


[n+(r’) - n (r*)] 

n=n ' . ® - - 

S 

we, get 


(5) 

(4) 


B[n^J - 1 ij ■" (5) 

How since the trial distribution corresponds to density 
of electron in the solute cell A after the size differoice has 
been compensated by an increase in energy and n-g in solvent cell 
B then E[n^]- - Eg corresponds to the negative of electioohe- 
mical contribution to the heat of solution because the change 
from n^ to ng will correspond a charge redistribution for 
equalising the chemical potential thicughout the systan lead- 
ing to a charge transfer giving rise to a dipole layer in the 
boundary of the cell A. The charge transfer Q can be estimated 


as 
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Q = / - n.(5)Jdr = - / [n-Cr) - n+(r)Jdr (6) 

cell A - cell B ^ - 

So, eqn, (5) oan be rewritten as 

Eg - E[n^] -i Q(m-^ - lig) (7) 

Since Q oan be shown to be proportional to - ji^) so the 
attractive electrochemical contribution to the heat of solu- 
tion is given by 

‘^e.o = - (^‘A- 

Identifying Pauling's rule for heat of solution with eqn. (8) one 

can arrive at n as a measure of electronegativity. In principle 

\x can be calculated from a Imowledge of B„[n(r)]* 'Haile 

2 

calculating and \i^ Hodges Stott took a uniform density 
model of pure components and were tuaable to get rid of r- 
dependence of n, So they averaged this ^i^(r) over A cell to get 
\ij^. This is arbitrary and so in our study we have taken as 
the value of at the boundary of the uniform distribution 

where the phenomenon of charge transfer takes place. The 
ezpression of eq.n, (8) also has been ezanined in this s'tudy to 
judge the validity of Taylorseries expansion upto 2nd order 
used in arriving at it, 

Q 

Miedema and co-workers , in theirestimate of this (xmtri- 
bution are consistent with Pauling's rule but surprisingly 
they have identified electron nativity with the work function. 
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3 

In view or the result of Hod-jes and Stott this again is 
il.u^'p.riaing because work function is the sum of the internal 
chemical potential p and the potential due to the metal sur- 
face dipole layer , So, unless it is shown that this potential 
due to the dipole layer is proportional to the internal chemical 
potential their measure of electronegativity is inconsistent 
with that of Hodges and Stott . 

In 1966 Hume fiothery^^ revisited this area and came out 
with Some more interesting details: (1). A tendency was observed 
for atoms of elements near the ends of the short peiiods and 
B-subgroups tocomplete their octet of elections, and a 
similar but less marked tendency to fill the d-shell in the 
later transition metals (2) Atoms giving rise to almost pure 
p— bonding may not enter structui^s involving other type 
bonds, A .high value of difference in electronegativity may 
overcome this restriction. The last observation cannot be 
commented upon in this study because of difference in the 
approach, 

1,1.3 Density Mismatch effect ; 

This factor is an artifact of the theoretical invest!- 
gations of alloys by Hodges and Stott"^. When the uniform den- 
sity of electrons in solute cell A that in solvent 

cell B (n^) readjusts in such a way that there is no further 
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charge transfer, at the houndary between these two cells a 

density step exists. This , due to the gradient tern in kinetic, 

exchange and correlation energjjwill give rise to an unphysical 

singularity in energy. So these authors have taken a small 

localized density component (u ) of exponential nature over 

and above the uniform distribution from both sides of this 

discontinuity. From the expression of chemical potential, with 

equation 

the help of poisson’s ^ . they find out the exponents of n^ 
and estimate the energy change due to this n^. 

This procedure is suspect atleast on two grounds (1) The 
component n is notdet emined variationally from the minimal 
property of the ground state energy (2) The reference state of 
no charge transfer with a step in uniform distribution is mis- 
leading because these authors get this state with an averaged 
chemical potential where its r-variation is smeared out. This 
may prove unreliable in defining a no charge transfer state. 

In this study we have not been able to separate this 
term in the strict sense. The comparison of the step model and 
the continuous model gives a rough idea only because in all the 
cases we found that the continuous models arrive at a minimum 
energy configuration by not only relaxing the step of the step 
model but producing a shift of the step also, hgain the value 
of the slope of the density distribution (-v ,7 to 3 ) B)»all 
indicating a significant charge transfer and consequently an 
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increase in the dipole energy. 

1.2 Eigid Band Model 

It assiinies that the eneigy surfaces and the density of 
states curve of the solvent netal renain unchanged on alloying, 
the only effect of the addition of the solute netal being, 
if its valencyis greater than that of solvent, to add ele- 
ctrons to the solvent band resulting in swelling of the 
Bemi-surface (P.S.) and filling the density of states curve 
to higher energy. 

By assuning the P.S. of- noble metals nearly spherical 
this theory provided an explanation for Hune-fiothery^^ 
observation of the solubility limit for f.c.c. a-phase at 
electron atom ratio of 1.36 in the Gu— Zn systen which 
coincided with the peak of the density of, states due to P.S* 
touching jhe zone boundary. Similar explanrtions were provided 
for the fixed electron-atom ratio for the limit of p and 
Y -phases in these alloys. But Pippard's eaperimental obser- 
vation negated this explanation by showing that P.S. has 
already crossed the 1st B.Z. after filling it up. In an 
.attempt to revive this explanation Hume-Eothery and Eoafe^ 
suggested that the limit of a-phase boundary is associated 
with the second peak of density of states where P.S. is 
making contact with the square faces of the zone. The 
defenders of this explanation ha^e further pointed out 
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inconsistency of Pippard's result with the electronic specific 

heat data which was interpreted again in tenos of the rigid 

band uodel. There has been attempts by several investigators 

12 

to clarify the situation from different points of view 
with limited success but still it remains an unsolved problem 
in alloy theory. 

In this investigation none of the above mentioned as- 
pects of rigid band model will be examined. This brief note 
originates for the sake of completeness only. One other 

experiment need be mentioned in this context - that of soft 

13 

x-ray spectra in Mg-AL alloy by G.A. Sooke . The spectra 
of Mg and A1 atoms in the alloy showed different band widths 
very similar to pure Mg and pure Al. This stands in glaring 
contrast to theassumptions of this simple rigid band approach, 

1,3 Priedal’s Model^ 

In a solid solution the perturbing potential produced 
in a gas of valence electrons by the solute atom is deter- 
mined in this model by solving Poisson's e(in. in the 
Thomas-Permi approximation, 

T = AS (9) 

p r 

where q* = 4 iiNq (Ep) and N^(^) is itLe density of states 
at the Permi eneigy Bp, is the valency difference between 
the solute and the solvent,, and r is the distance from the 



15 


solute. This expression of shows polarisation of the 
electron gas to screen the solute potential making it loca- 
lized. Since the solvent away from the solute atom is 
unaffected, the I'emi level ®p remains the same unless the 
impurity concentration is so high as to cause mutual 
interaction. The change in the density of states due to pola- 
risation has heen studied in great detail . 

The rigid band relation is found to be valid when is 
spherically syiametric provided it varies little within the 
wave length of the electrons wildi respect to the energy B-7 . 

Ir 

This is true when 2-7^ is not too near to the limit of the band. 
This approximation is also quite good at large r.. 

It was also shown that for strong localized perturbation 
bound states may occur between the energy bands and are sub- 
tracted from these bands. As these bound sta'tes are formed 
energy levels move out of the bottom of the conduction band 
with the topof the band moving down towards the Peimi level. 

Inspite of this detailed study of the perturbed electron 
states and th.eir polarisation, the absolute value of the 
heat of solution for dilute alloys could not be found out 
within this scheme. As a result this study finds little 
use in this investigation apart from one vital conclusion of 
Priedel^ to suggest that the wave functions will always try 
to adjust themselves in the alloy such that the solute atoms 
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are very mu oh like those in pure solute metal. Heine and 
14 

Vfeaire look at it as the quantum analog of the classical 
theoron that the energy density of black body radiation 
per unit frequency range per unit volume is independent of 
the boundary conditions likesize, shape or nature of the 
container. This is in complete consonance with the results 
of softs- ray spectra of Mg-Al alloys of Booke^^ as pointed 
out in the last section and it forms the starting concept of 
the approach' V taken in this investigation i»e, the solute 
atom goes in the alloy as a piece of metal of pure solute. 

1.4 Pseudopotential Perturbation Theory 

Since atoms as they exist in pure metal is a very 
relevant concept in alloy so we start with a discussion on 
how the pure metal is viewed in this formalism. The densi-ty 
of electrons in a pure metal is conceived to consist of an 
uniform electron gas with a weak periodic component superim- 
posed on it due to the p er'tTurbation produced by the pseudo- 
potential of the ions embedded in this uniform electron gas.. 
The energy of the uniform electron gas under pseudopoten— 
tial is the structure independent term and accounts for al- 
most ninety percent of the cohesive energy in many systems 
indicating that the structure dependent term resulting from 
periodic component is quite weak "bo be justifiably treated 
by per'turbation theory. But we find that “the metals with 
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higher valenqy and lower core radius (or ionic radius) the 
structure dependent term increases multifold in this linear 
response formalism, Dagens, Basolt and Taylor^^ have calcu- 
lated the charge densities and the interatomic potentials 
for Li"^, Mg and Oa'*"*’ ions embedded in uniform 
electron gas of corresponding metallic density self-consistently 
and found that the non-linear contributions are quite impor- 
tant. These authors tried to simulate this contribution 
within the firework of linear response theory by adjusting 
parameters of model potential but failed, Tet the success 
of this formalism in the calculation of different properties 
of pure metals stands out as a clear justification of its 
applicability. But in the study of the vacnacy or alloying 
in metals the perturbation may be quite strong even resulting 
in the formation of quasi-bound state and hence the second 
order perturbation theory may not provide a satisfactory frame- 
work for the study 'of the alloys. 

The alloy is conceived as independent atomic cells of 
solute and solvent . uniform electron gas of corresponding 

metallic densities. These cells are coupled by periodic 
electron densities of the alloy duly modified by the pre- 
sence of solute and distinct in that respect from that of 
solvent. The picture of coupling becomes physically trans- 
parent if represented in terms of interatomic potentials. 
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Here we present a brief accoimt of this formalism following 
Inglesfield*^'^^ mainly. 


1.4.1 Alloy potential : 


If^G-, conoentration of the solute with pseudo potential and 
concentration of B, the solvent with pseudo potential 
Vg is made into an alloy then the potential of the alloy 
system can be written as 


V(r) 



t I 

i(B) 




(10) 


where r is the space co-ordinate and represents the co- 
ordinate of the ions, i(A) indicates the sites occupied by 
solute A atoms and i(B) are those of solvent B atoms. 
Taking Pourier transform of V(r) we get 


V. 

ion 



1 I 

“ H i(A) 


-iq.E^ 

e v^( q) 


-Hi ^ 

^ H i(B) 


-iq.iL 

e ^ VgCq) 


= S(q) V -H [0 - (I-O)S^] (v^(q) - v^(q)] (ll) 


A-* ^''B' 


= i I 


where 3 (q) = ^ 4 e is the stiucture factor and 

, T y -iq.H. 

“ 1 i(A) ® H i(B) ® 


( 12 ) 


are the partial structure factors. The potential is now 
written in teimis of a mean potential v = 0 v^ -h (1 - 0) v^ 
and an alloying potential v^ “ ^A “ ^B* alloy 

can be seen as a periodic array of mean potential in the 
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uniform electron gas intermized with, a scattering alloying 
potential, 

1,4.2 Total eneigy : 

The energy of the uniform electron gas or the structure 
independent energy per atom of alloy is given hy 

V = C +Ug (1 - a) (13) 

where and Ug are respectively the energy of the unifoimi 
electron gas of densities corresponding to pure A and pure B 
respectively. 

Now the hand structure energy per atom of alloy evalua- 
ted by second 'order perturbation theory is given by 


= I 


'BS 


(llki*-Ik+ql 


(14) 


where V(q.) = V^^^(q)/e(q) 


gives the screened potential with s(q) as the dielectric 
function. To get correct value of we have to subtract 
the electrostatic self energy of the electrons ;diich have 
been counted twice in the above e3i)ression and finally 
we get 


Ub3 = I V(q) ^ e(q) x(<l) 


(15) 


where 
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where x-(<l) = | k< ^'^2 


Now using expression (11) in eqn, (15) and writting the 
in terms 

different terms,/of their Fourier transforms we get- the total 


hand structure energy as 


^BS - 


+ 0 * 


(BB) (AB) 


(16) 

( 2 ~t2k ) ^ 

where i2(j,(F) - J — + — ~ i [v(q)^e(q)x:(q) Je^'^'^d^g (17) 

is the mean atomic potential and 

(4)*E(<i)i(4)]e^1*®d3|l (18) 

4^ 

is the alloying atomic potential, and are valencies of 

A and B respectively. 


There will he an additional electrostatic term for assuming 

podvhedra 

spherical cells in place of Wigner-Seitz^Pg Pe^ atom of the 
alloy 


^ “a ■’* Rg 


(19) 


where R^ and R^ are atomic radii of ,A and B respectively. 

These oc^ and oc^ depends on the structure of A and B and is of 
the order of ,004. Now the total energy is given hy 


B = n +ffBa+''B 


(ao) 
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and tiie heat of solution per atom of alloy at 0 K is 
obtained by subtracting the energy of respective amounts of pure 
A and B alloyed and it can be expressed as 

AH = A U + A Ugg + ^ Ujj 
where, AU = U - OU^ “ (l-C!)Ug 

‘'^BS = '^BS - 

iUg ^ TJj, - 0l4 - (I-O)U^ 

Apparently AIJ and AUg become- zero. But if we write U 
and Ug not in terms of pure -.metal atomic radii but in terms 
of alloy atomic radii then these terms are nonvanishing. lU 
the three terns of eqn, (20) can be written as a function of 
alloy Tg, i.e. the radius of the volume occupied by an elect^n, 
then AH can be variationally determined with respect to r^. 

Prom this rg, alloy lattice parameters can be determined to 

check Vegard’s law. 

1.4.5 Interatomic potential '• 

in Seotion 1.3.2, eqn. (16) ve have expressed the hand 
struotnre oontrihution to the energy of -the alloy In terns 
of interatomic potential. So it is in order to write a 
short note on it mentioning the uncertainties and limitations 


involved. 
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The second term in eq.n, (16) involves an integration in 

the reciprocal space with the dielectric function appearing 

along with other teims in the integrand, Johnson^^ reported 

17 

that with Harrison local pseudo potential for Aluminium 
he failed to produce any minimum in the interatomic poten- 
tial, in the vicinity of the first nearest neighbor by 
using Hartree dielectric function Harrison confirmed this 
and contends that the minimum appears if the non-local fiill 
ps eudo— pot ential is used, Ashcroft claimed that the 
interatomic potentials are quite sensitive to iincertainties in 
exchange and correlation and to dielectric function e in 
the region 2 < k < few kp. It is possible in some instan- 

ces to eliminate the principal minimum in interatomic poten- 
tial by ignoring exchange and correlation. The effect of 
density dependence of the dielectric function does have a 
Substantial effect on the inteiatomic potential. The calcu- 
lation of Dagens et al^,^ shows marked effect of nonlinearity 
of the response function in the first nearest neighbor 
region. 

From the above it can be clearly inferred that these 
calculations are extremely sensitive and so their reliabi- 
lity and consistency cannot be commented upon unless a clear 
understanding evolves regarding the factors involved. 
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1* 4»4 Ohoice of model pseudo::'ot ential ; 

In the calculation of interatomic potential in alloys 
the choice of the model pseudopot entials for the constituents 
is the starting point. From the scattering approach this 
model potential can be viewed as one which simulates the 
scattering of the band electrons correctly and hence it is non- 
unique* These models, in general, contains free parameters 
to be adjusted by comparing with some observed properties. The 
Abarenkov, Animalu and Heine (AAH) model is fitted with the 
spectroscopically observed atomic energies whereas the 
Ashcroft empty core model is fitted with the resistivity data 
or the lattice parameter data. These are the two -models of 
relevance to this study and so we briefly describe them. 

From Herring's OPW method of band structure calculation one 
can write 

. v(OPW) = ^real^^^^, ^^^c^ ’*'c 

c 

where c indicates summation over the core states with the 
energies and the wave functions This equation means 

that orthogonality of the valence states to the core states 
can be taken care of by an equivalent transformation of the 
potential as expressed in eqn. (22), 
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Sincf^ v^g^(r) is negative and B^'s are negative, the two 
terms of eq,n«(2S) have opposite sign and cancel to an extent 
inside a radius where ip-(r) is substantial. This 

radius is called the core radius inside which the core-electrons 
are confined. The second term of this eqn, (22) is nonlocal 
in nature and if we write it in the form 


v(r) = f(r) + I f3_(r) 


(23) 


where is the projection operator separating out the 
corresponding 1-component from the wave function. Prom the 


identity I P = 1 we write 
1 ^ 

0 = f - f 
av av 




(24) 


Adding eqns. (23) and (24) we get, 


v(r) = [f(r)+f3^] + I (f;j_(r) - 


(25) 


If f is large and positive like the f, then both the Bracks- 
eted ternB become, small. So even if the cancellation within 
the core is incomplete from eq,n, (22) a shift of origin by 
proper (Choice of may make the non— local term in the 
potential negligible in most of the cases for simple metals. 
Prom this justification Ashcroft proposed a local empty 
core model for pseudo— potential of the form 



v(r) = 0 for r < r 

G 

= -Z/r for r > r 

O 
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(26) 


Here r is the parameter to be fitted with an observable 
property of the system under consideration. In our study 
we have used this model for our parametric variational 
calculation of the heat of solution of the alloy. In our 
pseudo atom calculation reported in the last chapter. We 
have used (for reasons stated there) a restricted form of 
AAH-potential optimised by sha^'^to make it smoother in real 
space and cut-off faster in reciprocal space. This is origi- 
nally a non-local potential but we have used a local form 
by taking only the 1=0 component. So, the form be<>omes 


v(r) = -A„ for r ^ r. 

= -Z/r for r > r^ (27) 

Shaw^^ has proposed a value of equal to -Z/r^ which 

flgflT n makes it one parameter potential if higher 1 ■ non-local 
contributions are neglected. It is eqn, (27) along with 

— — 2/r which we have used in the last part of our study. 
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Density functional Method in an Alloy 

This method expresses the ground state energy of a many 
electron system as a functional of its correct charge density* 
This charge density in turn can he found out hy using the 
minimal property of the ground state energy and the functio- 
nal relation. The formalism is exact in principle although 
it involves some appixDximations due to our lack of knowledge 
of the exact functional relation for the kinetic, exchange and 
the correlation energy. This chapter first describes this 
formalism and then approximations normally used. 

Then we show how the density functional formalism can he 
cast as a neutral pseudoatom theory. This has heen done hy 
Dagens^ for the pure metal case starting with some arhitraiy 
parameters. We use two types of formulations (i) Parametric 
variational where the charge density is parametrised and 
the parameters are evaluated hy minimising the energy, i.e. 
to get the ground state energy. This scheme is described in 
detail in the next chapter, (ii) The Eohn-Sham formulation 
where we get Schrodinger like pseudo one electron equation 
from the variational principle. This scheme will he dealt 
with a view to apply it to pure metals and alloys in Chapter 
5. Here we deal with the basic formalism in brief and 
discuss how it is applied to the case of alloys in general. 
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The idea of applying this method to alloy came from its 
successful application to the problem of metal surface where 
the situation is q,uite similar* In the metal surface the 
electron density changes from that of the bulk to zero ower 
the surface layer whereas in the alloy it changes roughly 
f 2 ?om the bulk value of the solute to that of the solvent 
over a similar distance* For the metal surface both 'the 
parametric variational method (in which calculations were 
carried throxigh first by Smith ) and the Ebhn Sham method 
(worked out by Lang and Ebhn”) give good results for the sur- 
face energy, work function etc. We thus ezpect, a forteriori, 
that the methods should be quantitatively successful for the 
alloy p 2 ?oblem also. 

2.1 Density Functional Fomalism 

In 'jhis section we will review the general aspects !#f 
this formalism leaving out the particular application to 
metals and alloys which forms the bulk of this study. 

Hohenberg and Eohn^^ (HK) proved a remarkable theorem 
which states that the ground state wave function of a many 
particle system, if non-degenerate, is a unique functional 
of the particle density* This theorem implies the exi- 
stence of a universal ground state energy functional of the 
external potential and particle densily. The eneigy 
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fimotional is minimum for the true particle density, so 
that the true ground state energy and particle density may 
he determined by application of a variational principle 
involving the pa2rticle density alone. The original proof was 
for a local potential and was extended for non-local poten- 
tials by Gilbert^^. 

In order to apply the theorem, one must be able to 
construct the universal functional. It is accurately known 
for the limiting case of a unifom high densi-ty electron 
gas An approximate functional form for slowly varying 
electixin gas^^*^^ is known even for metallic densities^^. But 
serious questions have been raised regarding the validity of 
the inhomogeneous term in view of the fact that the density 
gradient expansion near a nucleus where the density gradient 
is very high, can be shown to be divergent series and so need 
be Summed to infinite order^^. Local spin density theories 

26 27 

based on HE! theorem have been applied to atoms and molecules 
inspite of these troublesome questions with surprisingly 
good results for many applications. In all real molecules 
and solids the external potential results feom one or more 
nudii which produces a very inhomogeneous electron density 
distribution. But in solids if one is interested in 
problem concerning the valence electmns which have a more 
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nearly unifoxm density, especially if one uses a pseudo- 
potential to smooth, out the density in the core region, this 
formalism provides a tool of enormous strength. Our appli- 
cation lies in this domain, 

2,1,1 HE Theorem t 

Let us ooneider an electron gas under the influence of 
an extemal potential v(r) and the mutual ooulomlj repulsion. 


The Hamiltonian has the fom 

H = T + V + U (1) 

where (in atomic unite) 

T = i (v) dr (2) 

Y - ! v(r) 4'^ (r) it»(r) dr (3) 

U = — 1 ? (r) (r') (r’ ) t); (r) dr dr’ (4) 

M 

The electronic densily of ground state is given ty 

n(r) “ ( J f ’I' (r) (r) $ ) (5) 

which is clearly a functional of v(r). 


How it has to he shown that conversely v(r) is a 
uniq.ue functional of n(r) apart from a tiuvial additive 


constant. 
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Proof 1 j-ssuDie another potential v(r) with ground state f 
giving rise to the same density n(r),ii $ since they 
satisfy different Schrodinger eqn. unless v'(r) - v(r) = 
constant. Hence if we denote the Hamiltonian and the ground 
state energies associated with $ and ^ ' hy H,H' and 
respectively, from the minimal property of ground state we 
get, 

B' = (J',H' ) < ($ ,H' $ ) = ( i ,^+ V - V) i>) 

or, B' < B + J [v’(r) - v(r)] n(r) dr (6) 

Similar arguments by interchanging primed and unprimed 
q.uantitioB will lead to 

E < B' + / [v(r) - v(r)] n(r) dr (7) 

Addition of eqns« (25) and (26) yields 

B - B' < B - B’ (8) 

which is an absurdity. Thus v(r) is (to wit h i n a constant) 
a unique functional of n(r), since in turn v(r) fixes H we 
see that the full many particle ground state is a imique 
functional of n(r). 

Since $ is a functional of n(r), so are evidently 
the Iclnetio and interaction energies# So one can define 


a functional 



( 9 ) 


i'[n(r)] = ( , (T-rtJ)i ) 

■wiiicli is 8, univ 6 2 ?s al functionajL ^ valid for ary mmlsr of 
particles which is again a functional of n(r). 

For a given potential v(r), the energy functional is 
defined aB 

= / v(r) n(r) dr + F[n] (10) 

For the correct n(r), l^[n] equals the ground state energy. 

Now it has to he shown that F^[nJ is a minimum for the 
correct n(r) thereby yielding the ground state energy 
through a variational scheme. The admissible class of fun- 
ctions are- restricted by the condition 

N[n] = / n(r) dr = N (ll) 

For a system of N-particles the energy functional of ii> ' 
is given by 

= (J’fV $*) + $') (12) 

It is well known that e^[ has a minimum at ‘the correct 
ground state J relative to arbitrary variation \mder the 
condition that the number of particles is constant. 

For a given external potential v(r) 
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“btiL help of eq.ns^ (9) xnd (12) one can. write 

/v(r) ]a'(r) dr +F[n'] > / v(r) n(r) dr + P[n] (15) 

where n'(r) corresponds to state , Equation (15) esta- 
blishes the minimal property of the energy functional with 
respect to the density. This provides the basis for the 
variational scheme in this formalism, 

2,1,2 Local density appicximation i 

The basic idea is to assume local uniformity in an 
infinitesimal volume element of the co-ordinate space and 
use the well known expressions of linetic energy, exchange 
and correlation energy for the uniform electron gas within 
that volume element added with the first term in the gradient 
expansion of the energy functional. So it can be under^ 
stood t**at such approximation might work well only when the 
density is slowly varying in space. 

Prom eqn, (10) we can further separate out the classi- 
cal coulomb energy and write 

E_[nJ = / v(r) n(r) dr + i f f dr dr’ + &[n] (14) 

~ - |r - r’f 

where G[nJ is another universal functional containing kinetic 
exchange and correlation energy. 
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How ©[nj can b© ©zpandsd in "Geniis of a gsnex'Eil gnadieni 
©zpansionj "blie firsi "two "fcomiis of wbioh. we genorally us© 

G-Cii] = I S^ir) dr + / g 2 (r) |^n(r) dr + (15) 

In this approximation the first term can be written as 


! S^ir) dr = / t(n) n(r) dr + / e^^(n) n(r) dr 
where t(n) = ^ 


( 16 ) 


r (n)‘ 


= s^(n) + e^(n) 

2J.QV / z' o' ' 


e,(.) = - I = - 0^ 




For £q(ii) different expressions are available starting from 
28 

that of Wigner ° to the most recent one by Gnnnarson, Lund- 

2Q 

qvist and lundivist . For reasons to be stated later, we 

30 

used the ©spiressions provided by Hozieres and Pines-^ . 


Sq = -0.0575 + 0.0155 In rg(n) (17) 

One can identify t(n), IcLnetio 

exchange and the correlation energy per electron of a uni- 
form electron gas of density n. For further convenience it 
has also been expressed in terms of local electron liquid 

parameter re,(n), i.e. the radius corresponding to the 

s 
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sph.6r© Gon'taiiiiiig on© ^ecbron. In oui* step model we use 
this approximation for G without the gradient tem to he 
discussed in th© next section, 

2,1,3 The Gradient Term : 

Th© contribution to the gradient teim can be separated 
into two parts (a) that coming from kinetic energy and (b) 
that coming from exchange and correlation energy, 

82 ^^') = [n(r)] + (18) 

8 n(r) 

X is a paz*amet©r in the kinetic energy contribution to the 

eq.n, (18) which was originally proposed by Von Weizsacker^ 

with X = 1, This form has been extensively used in th© 

32 

literature. Jones and Young have compared the response 

gradient espansion 

fuhotion obtained from the truncated/ with variable co- 
efficient X to th© Lindhard response function. They found 
that choosing X = 1 yields the correct response for pertur- 
bations whose wavelength of importance is small with respect 
to the Fermi wavelength, while x = ^ is appropriate for 
perturbations whose wavelength is large with respect to the 
Fermi wavelength. Thus, choosing x = ^ we woiild expect 
the induced screening charge to be relatively accurate at 
distances far from nuclei or ion but relatively inaccurate 
at distances near to it« HodgeS et al have recently 
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dstermiin. i tJie nexrt term in the gradient ex^tansion eq.n«(l5) 
as contribution from kinetic eneigy, 

[(2Lfi)2 - I (?^) (^)2 + i(9a)td? 

(19)’ 


But unfortunately it was found by Wang^'^ et al that inclu- 
sion of this term increases the deviation to +5 percent from 
kinetic energy estimated by Hartree Pock Method for noble 
gases. If this teim is dropped then it remains within less 
than il.O percent. So in this study we ignore this term. 


The exchange and correlation contribution to the gradient 
term has undergone an intensive investigation in recent years 
through the studies of Sham, Herman, iiasolt et al and 
Rajagopal and this author^^. The results of the last two 
investigations are the latest and agree in general with 
each other. Since the results of this investigation have 
alreaciy been published we will include the details of this 
work only in appendix. In this study an . exact and a varia- 
tional solution of the electron scattering yortox e<iuation 

has been obtained. Prom these, one can arrive at the 

( 2 ) 

static density correlation function to find out g^^ , i,e, 
the exchange and correlation contribution to g 2 of the 
gradient term of eqn, (9). apply this general expre- 
ssion of g^^^ to Yukawa type of interaction of electrons 
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with a screening parameter^. In the large^limit the varia- 
tional result coincides with the exact result but for 
5 ■*■0* i*6* strictly coulomb limit there is a logarithmic 

divergence for the exact resiat. In place of using Thomas- 
Permi expression for E if one uses a self consistent 5 
from the longitudinal dielectric function our results agree 
with that of Sham in the high density limit and that of 


Kleiuxan for iskj, o (low density). 


« 


( 2 ) 

xo 


C (rg) n 


- 4/3 


A plot of o(r-) vs r- has been provided by Basolt^^, In our 


S' 


‘S 


continuous density model of alloy we use the local density 
approximation along with the gradient terms discussed above i 


2,1,4 jSxaot Kinetic Energy Pormulation t 

In this formulation of the density functional method 
Kbhn and Sham^^^found a way to avoid the local density appro- 
ximation for the kinetic energy as folloijs : 

Suppose for a given v(r), ii(r) is the correct density 
corresponding the ground state eneigy from e<in, 

(14), If we allow an arbitrary variation of n(t’) subject 
to the condition that / 6 n(r) dr » 0 the minimal property 
of the ground state energy demandi that 6 E^[n] = 0 or in 
another way it can be expressed as 
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r 

J 6 ]i(r) , [ _ 0 ( 20 ) 

6n(T‘) 

for arbitrary 6 a(r) subject to the restrictive condition 
stated above* c is a Lagrange multiplier# 


6 \ 

6 n(r) 

where 

and 


« ♦ (r) 
♦ (r) « 
(n) - 


a TgEnJ 

T(r) + / dr’ 

^ (GEn] )) = 

6n(^) 



( 21 ) 

( 22 ) 

fe (n.„(n)) 

( 23 ) 


B_Q[ixj le the exoh£al£e and correlation contribution to the 
(J[nj and can be written approximately as 

with e^^[n] the exchange and correlation energy per electron 
of a uniform electron gas of density n. This gives a ther- 
modynamic meaning to in the sense that it is the exchange 
and correlation contribution to the chemical potential* 


Equation (20) written in operator form loots like an 
one patrtiole 3ohrodii>€«^ equation* 

[- i A* + C ♦(!') + |AjQ[n(r)] ] ♦^(r) » (^25) 

n(r) « J |^.(r)P 

i*a . I 


with 


(26) 
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where N i3 the number of electrone, ihe set of equations 
(22), ( 23 ), ( 25 ) and (26) can be solved self consistently 
to get e^'s and n(r) which in turn will yield the ground 
state energy as given below, 

B - ^ - -i- / / + j n('?)[e^^(n(r)] 

- ^ 3 ^oCii('r)]j dr (27) 
where the second term is due to double counting of electron- 
electron repulsion term. 

The above procedure has been used in this study in the 
last chapter to find out the cohesive energies of pure 
metals and the heat of solution of alloys from the point of 
view of a pseudo-atom approach which will be discussed in 
detail there, 

2,2 Density Functional Scheme and the Neutral Pseudo-a'tom 

Oonoept 

We briefly outline here the various ways in which the 
density functional scheme can be implemented for oaloolating 
the ground state energy of pure metals and eHoys* Many of 
these are best understood, in terms of a neutral pseudo— 
atom idea. We first discuss the pure metal case, and then 
the alloy. Details are given in <3hapters 5 and 4, 
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a) Pure metail : 

The most direct method, followed by Jauak et al^^ is to 
start with an external potential v{r) (periodic), solve the 
Kohn Sham one electron equation using one of the methods 
standard in band theoiy, determine energy levels and 

wave functions construct nC?) (=s 

and thence r^^^ (n(r)) and iterate to self-consistency* 

This has been done, with v(r). l — where 2 is the 

I ^ ^ ^ 

nuclear charge, by Janak et al for transition metals 

with great success* Olearly the method avoids all assumptions 
about the ionic potential: the only approximations made are in 
the density functional form. However, the method is oomputa- 
tionsilly prohibitive, and physically not very illxaminating. 

In general, the density n(r) can be written 

n(r) x: n + f 6n(r-H.) (28) 

0 j , J 

where n is the average density, and in a pure metal 
0 

/ 6n(r) dr = 0 (29) 

This form suggests the following two approaches. 

' 

i) Since u(r) is periodic 

iK.r 

n(r) = n^ + 1 Ug e 
£ 


(30) 
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S beiii« the recipi?ocal lattice vectors of the lattice. 

Thu-S in a crystal | the energy is a function of the discrete 
set of variables (n^,[n^]). One knows the functional foim 
for = 0 (unifom electron gas), and the form for 

small [nj] (power series in nj.). So, if nj/n^^ ^ 1, one can 
proceed as follows i 

ilhe H-K eq.n, (10) gives the ground state energy fun- 
ctional as 


B * f v(r) n(r') dr + P[n(r)] 


a / v(r) dr + I nf + P(n,n£) 

£ 


(31) 


where v('r) is the sum of the ionic potentials and S’ is the 
energy of interacting electron gas. If ng*s are small then 
one can ezi)and B in a Taylor series as 

P(n^,nf) E(n^) + i I n.J (• * 


K 


^•-3 nK 


)(ng= 0 


+ i i Ilf' n.f + 

£ 


(32) 


The term linear in n^ is absent because n^ is the mean valence 
density and so 6n(r) integrates to zero. 

Minimizing B with respect to n^'-froia eqn. (3) after substi- 
tuting eqn, (4) for B we get. 


n^; = - 




( 35 ) 
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We note t:.at n£^ is small if ^i) small znand (ii) the 

electron gas reflponds weakly. Putting this n^ from eqn, 

(5) into eqn* (3) we get the ground state eneigy 



J v(r) n^ dr + PCn^) - i 



(34) 


In linear response theory the last teim of ein, (34) is 
identical with the second order perturbation term in energy 
in pseudopotential perturbation theory. 


The above derivation shows the unity in two apparently 
diverse approadh and the separability of energy contribution 
due to the uniform and the periodic component in the domain 
of validity of perturbative treatment i,e, n^^' small. 


ii) The neutral pseudo-atom appmaoh. The idea is to 
develop an expansion for 6n(r-E.j), Clearly, the major role 
in determ. n.ing 6n(r-Rj) will be played by the potential in 
the Wigner Seitz Cell j, Since for r within this oell, the 
potential is strong. The effect of other c^s is veiy small 
because the i^rcic potential is nearly completely screened by 
the electrons I the residual effoct is weak and of short 
range. So, one determines 6nCr-Rj) non perturbatively in 
so far as the effect of v(r-£J) and the screening charge in 
the cell j are concerned. Such a programme can be implemented 
in two ways. 
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ii oc) Parametric method : 

We parametrize 6n(lp) suitably, and use a known density 
functional form, e*g* the local form plus gradient tenns. 

This method is non perturbative in the sense that 6nCr) is 
not a priori assvimed small j it is variationally deteimined. 
The method is discussed in detail in Chapter 3 and is applied 
there to some simple metals. It turns out that the effect of 
other cells on 6nCr-E^) is negligible, occuring only via the 
electrostatic interaction teim in E[ n J* 

ii p) Kohn Sham method i 


As stated above, one can obtain an approximate expression 

for 6n(r-R. ) by neglecting the effect of other cells. The 
J 

effect of the cell J (i.e, v(r-R. ) and the effect of electron 
density n^ inside cell and of 6n(r“R^")) is considered 
• exactly S i«e* we solve the Kohn Sham equation for the 

, obtain wave functions 

struct 6n(a>-R^) and iterate to self consistency. (See 
Chapter 5 for details). Prom the n(r) so obtained SX n" j 
can be calculated. Howevert here again, the effect c£ two 
centre terms is calculated perturbatively, i.e. 


appropriate V 


eff 


(r-R^) 



B[n. + '■ I 6n(^-R.)] 

0 j 0 , 


two center 


•i- ^ 6n(^R.) 6n(r-£) x n(r-Rj^) ^ 

^ ^ >^nCr-5^) f 


6n=0 


(35) 
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Sucsh. a ca^-culation was first carried tiirough by Da^enS, 

We present a systematic theory in Chapter 5 which avoids 
the introduction of arbitrary form factors as done by 

9 

Da^ens , and can be used to study the errors introduced in 
the approximations made. 

b) Alloy ; 

The problem of applying the density functional formalism 
clearly, in alloys, the procedure used by Janak et al for 
pure metals is inapplicable. However, if the ion potential 
(pseudo potential) is weak, the ground state energy can be 
calculated by perturbation theory as described in Chapter 1. 
The use of linear response theory to describe the effect of 
charge contrast — Zg 0) on the electron gas is questio- 
nable. However, since the effect of an impurity ion is 
strongly localized, we can u e the neutral i seudo-atom idea. 

Consider the case of alloy where neutral atom of solute 
A is dissolved in(the site at origin) solvent B. In the 
similar way as it was done for the pure metal case we can 
start with a rigid jelly of electron of mean valence density 
n-Q corresponding to that of B» Then we put in the origin the 
solute ion with the excess electrons or holes over the 
valency of B and the B ions at other sites. How if the 
electrons are allowed to relax the B— ions will be 
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independently screened but near the A-site two phenomena will 
take place (i) the screening of the solute ion A and 
(ii) the alloying i.e, attaining thermo-dynamic eg.uilibri\m 
by equalisation of chemical potential throughout the system* 

If n^ and 6 n^are "tbe mean valence density of electrons and 
the screening charge polarisation in pure A respectively and 
n^ and 6 n^ . ♦ are those corresponding to pure B then 

the charge polarisation of the electrons over n^ in the 
neutral A atom in the alloy will be different from ^A “ “b 
d n^ because of the effect of alloying. The charge distri- 
bution in the alloy can be written as 

n(r) = n^ + 6 n,^(r) + 'I- 6 n^Cr-R. ) (36) 

Note here that 

f 6 n^'!^ dr = - Zg in place of zero for the pure case. 

Zj, and are the valencies of the solute A and the solvent B 
respectively. In pure metals the charge polarisation due to 
screening gives rise to a weak interaction between the nearest 
neighbors ( ^ 4 percent of the energy of the neutral atom for 
Alkali Metals^ )but in the alloy it will not be small arouhd 
the solute atoms beca\ise charge transfer will cause 6 
to extend into the neighboring neutral atom cells wiiii 
magnitude depending on the difference in the chemical 
potential. 
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In £ parametric variaf.onal approach, r'ur attempt is 
to simulate the 6 n^-p^(r) by discontinuous and continuous 
functions in the step model and the tangent hyperbolic 
model respectively (See Chapter 3)* 

The etjuation (36) can be rewritten as 

n(T) = (37) 

where iig{r) = nj - n? © (Rp -|f r|) + ^ 6 nS(?-Bi) 

- 1^0 

= >^ 9 (Bb - I r||) + 6 

The atomic radius of pure B and its lattice vectors are 
given by Bg and R respectively and n^ is the density of 
elect2X)ns of the pure solvent as it takes pant in the alloy, 
i^^^lr) represents the electron densiiy in the solute atom 
as it exists in the alloy. Here we assume a spherical 
Wigner-Seitz cell for the solvent with a radius R^ in 
place of a polyhedra. This is quite a good approximation 
for close packed structure although the electrostatic eneigy 
estimated will be slightly different from that deteimined 
by the Bwald method for a polyhedra. This will be discussed 
in more deta i ls in Chapter 5. The advantage of this separation 

— 4 - 

of n(r) into ^^(r) and ^qq(^) that the heat of solution 
of the alloy can be written primarily as a difference 
between the energy of A-pseudo-atom as it exists in the 
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alloy and the energy of an A- pseudo-atom as it exists in the 
pure solid. Otherwise if the heat of solution is -written as 
the difference in the energy of i the constituent atoms to 
that of all the participating atoms in their pure solid 
state it poses a problem of evaluating two large energies 
and the errors involved may make the estimate of the heat of 
solution totally unreliable. 

In the Eohn Sham version of the density f-unctional 
scheme, 6n^(r) is obtained by self consistently solving 
on one centre Schrodinger like equation, with v(r) 
and the other coulomb potential terms arising from the 
electron charge 2^ composed of distributed in the WS cell 
with densi-ty n^ and 2^ - Zg distributed with densi-ty 
6 n^('r). Knowing 6 n^C?), one can calculate the ground 
state energy of the alloy (soe Ohapter 5? Section 2), 

2,3 The Heat of Solution of Alloys 

When we refer to the heat of solution A ^ this 
study it always refers to the partial heat of solution per 
atom (the equivalent of partial molal heat of solution in 
the conventional thermodynamic term). If n(r) is the 
charge density distribution in the alloy and B[n(r)J gives 
its ground state energy, the heat of solution per atom 
of Solute is given by 
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V = B[a(r) j - M (38) 

wliere M is tJae nuinber of solvent atons B with, valency Z-n 
alloyed with 1 aton of A with valency Z^, Eg and are 
the ground state energies per atoia of pure solid solvent B 
and Solute A respectively.' Ihe foim of eqn. (33) » although 
correct, is not convenient for calculating the heat of 
solution for the reasons mentioned at the end of the last 
section. So we tiy to write 

E[n(r)] = M Eg + (39) 

so that , 

C40) 

where may be thought of as the energy of the solute 

atom A in its configuration in the alloy. Now the problem 
reduces to partitioning E[n(r)] as given in eqn. (39). Then 
the evaluation of the heat of solution by eqn, (40) becomes 
quite straightforward. 


Following eqn. (14) one can write 

B[n(r)] = fv(r) n(r)dr + 1 dr dr’ + G[n(r)] + 

■ ■ ^ Ul) 


where is the electrostatic interaction energy between 
the ions and G is the universal functional of n(r') containing 
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"bhe kineti 3, exchange and "the correlation en 'ivgy functional 
as already defined in the first section. 

Here the potential in the v(r) is given by 

= v-^(r) +^1 ^ Vg(ia^.) (42) 

where is the potential of the A-ion at the origin and 

VgCr-H ) is the potential of the B-ion at the site with 
lattice vector • Although in this study we will use the 
model potentials for both v^ and v^, it is not a feature of 
this formalism by any means, For a angorous calculation one 
should use potentials constructed forom the fiiet principles 
computations of the atomic stiucture as is customary in the 
different band structure calculations. 


Following the discussion in the last section the charge 
density in the alloy can be written as 


n(r) = + ^oc^^^ 


(43) 


Putting equations (42) aJih (43) in the eqn, (41) one gets 

E[n(?)] = ' 


0 


G[nBJ + X 

■ -I 

dWo 


V 3 (r^R^) n^{r) dr + i :~7T i+-- 


l?-r' I 





49 


+ dr 




loc' 


■¥ -V 


I I 


dr dr’ + G[n(r)J - 


- sender)] 


+ / Tj^{r) HjCr) 4? + J" 




-^1 

r-r' 


dr d?' 


+ { v^(r-fi ) 

-1^0 •“ 


“lo 



+ Z 


'A 


1^0 



(44) 


Since is tiie correct ground state density of the 

assembly of M atoms of pure solvent, the first four teims 
in eq.n* (42) gives MEg* The nezt four* terms are primarily 
one centre contributions originating from the pseudo-atom A 
as it exists in the alloy. It also includes a small contri- 
bution from the effect of environment of B on the kinetic, 
exchange and the correlation energies of the pseudo-atom A. 


The last four terms represent the electrostatic intera- 

■ • 

ction between the pseudo-atom A with the surrounding B- 
pseudo-atoms and so are two centre contributions. Since 
these terms contain long range functions we rewrite them 
in a form in which the long range parts explicitly cancel* 
These are computationally suitable for an empty core 
pseudopotential both for A and B, 



dr 


(45) 
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^3^ ') ~ 0 by definibior vithi n "blie or re radius of A. 

Again, 


■' ^ ^ 


^ ^ M ^"D ( ^ ^ r 

dr dr' = dr' / . ni„„Cr) dr 

\IrM . 0 






1| r-r'il 

~ jl-t ^ J 

'' Ir'l . . * ^ 

„ » , C'^loot'^’) - -f - 

= ^ ^ °° ' dj 

( 46 ) 


dr 


4 ^ 

\ 




where Q^Cr) = [ dr and '^oQ(r') = / ) d? 


. 0 


oc' 


-oc' 


Adding eqns, (43) and ( 44 ) we get 


^ T^(r) n3(r) dr + f dr| = 


+ 3 


S' 


( 47 ) 


where we have neglected 6ng as defined in the last section and 
so we can write 


Qg(r) s= 0 for r < Eg 

= |2E (r^-E^) ng for r fig 


( 48 ) 


Tile justification for ignoring; 6n° or otherwise 

'*^^‘-‘-1 be 

provided by the results of our calculation of cohesive 
for pure elements assuming a cosine function for dh? 

th 

some parameters to be determined variationally from th 
minimal property of the ground state ener^ 

. XII e estiQat 

of the amplitude of this cosine function will show how 
we are justified in the above assumption in different 

The eqn. (47) is clearly in a computationally suitabi 
form as both the integrals reduce to zero beyond 
distance. 


energy 

Le 

^ A 

e 

ouch 
cases. 


The last two terms in eqn. (44) can also^recast as 


j. 




with 


(i) = -( 


loo 


where r® is the core radius of solvent B and v /'v.\ 

° oc'’^-' is a 


located function. The first integral in the l.H.S, of 
(47) can be written as 


e<lh. 


= - I > - aj 

Vo JTi ■ 


/d? 


^Ooo 


d) 


,{ ^3 I ©( t - t )>- 2 , 2 




1=0 




‘B 


"B 


(50) 


U.T. ^ 






i 


lx .5?IH 
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Tile last t- im of eq.rL. (50) oaniels with, the last term of 

eqn. (42) and finally taking eqns. (-41), (47), (49), (50) and 

(40) together we can write the final e 2 j)ression for the 

heat of Solution A as 

A 


dr^dr» + 

I r-r’{ 

+ G[n('r)] - G[n°] 


.0 I ar . / 


+ /n^^^(r) dr I 0((rl - ftRj!}}- 


(51) 


This expression contains all short range tems only, 
by its construction is short range. So, the first two terms, 
the sixth, seventh and the ninth terms of eqn. (51) is 
shoirt range. appearing in the fifth and the eighth 

terms goes to fast enough for distances' somewhat greater 
than Rg rendering both of these tems also short range. 
G[n(r)j - G[n^J if expanded also leads to residual shoirb 
range terms making the eqn. (51) totally short range. 
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We use this expression for the heat of solut_on for the 
parametric 'variational scheme described in the next chapter. 
But for the Kbhn-Sham scheme we use a different version of 
the same formulation without assuming Sn^ negligible to 
be described in Chapter 5. 



Oliapter 3 


iaranietrio Tariational Method 

Here "we describe the ground state density of the electrons 
in the alloy by a function involving free parameters in it and 
compute the ground state energy from its minimal property 
with respect to the density variation i,e, the variation of 
these free parameters. From this alloy ground state energy 
one Subtracts the ground state energy of the components to 
get the heat of solution of the alloy. 

It is obvious that ary approach to the heat of solution 
of alloys has to be preceded by a treatment of the ground 
state energy of the pune constituents within its scheme. 

Also, for understanding the physics of the alloy problem it 
constitutes logical path to follow. So we start this 
chapter with a disciission on the calculation of the ground 
state energy of pure metals - its pTOmises and pitfalls, 
to provide for the necessary orientation to go to the alloy 
problem. 

The discussion on alloy problem will start with a 
a uniform density model ^or the components with different 
valencies resulting into a step in the alloy electron 
density. The two parameters in the density will be the size 
and the position of the step one of which can be determined 
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from "tile ions ervaf ion of chaise* The ofher parameter acts 
as the free one with respect to which the energy is mini- 
mised* This is essentially a one parameter variational 
model for the alloy* This model is quite insufficient for 
the treatment of the alloy for reasons to be mentioned in 
this section but it is an extremely simple model bringing 
out some aspects of the alloying phenomenon seen from e35>e- 
rimental studies so far* 

In the above model the gradient teim in the density 
functional expansion of the eneigy will have to be neglected 
because the discontinuity in n(rj at the step will make 
this term singular. To take into account this term (which is 
by no means negligible) we model the density with a conti- 
nuous tangent hyperbolic function with three parameters one 
of which ill be detemined by charge conservation* Thus 
this is a model vith two free parameters with respect to which 
eneigy will be minimised* The distribution being physically 
quite similar to the step model, this can be thought of as a 
model which allows the elections to rela:x from the step 
configuration with the help of the driving force provided 
by the gradient term* The inadequacies of this model will 
be brought out in this section for a clear understanding of 
the directions for future effort* The discussion on the 
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conipu'ba'tioiial resul1;s of 1/1168 6 Eiod6ls will t)6 posi:poii6d 
fill the nezt chapter where we will discuss the systenatic 
trends in the heat of solution of the alloys, 

3»1 Oalculation of ground state energy of pure netals 

Prom e^n, (14) of Ohapter 2 one can write the ground 
state eneigy of a system with the ground state density of 
electrons n(r) as 

B [nj = /v(r) n(r) dr + i J / ^ dT dr* + G[n] (l) 

S . r - rM 

where Gr[n] is a universal functional of 'density n{r) containing 

the kinetic, exchange and correlation energy contributions. 

This equation shows that the problem of calculating is 

. S 

essentially one of arriving at a correct v(rj and n(r) for 

the ground state of the system, Por v(rj we choose the 

1-shcroft local empty core moael pseudopot enoial as it yields 

of 

the cohesive energy valu^ the real systems quite in agree- 
ment with their experimental value. This indicates that 
empty core pseudo potential may provide a sufficiently good 
model potential for our purpose, 

Uow the next task is to choose a correct n(r) with 
Some variational parameter to be determined by the minimal 
property of ground state eneigy. Our choices for this 
purpose and their justifications are described below. 
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energy hae been written as a function of atomic radius which 

is related to the size of n^ for a given Z. The expression 

of B when minimised with respect to B gives the ground state 
s 

energy. By comparing this value corresponding to 

minimum and the B with the experimental value we can get 

s 

an idea as to how good is this uniform density assimption 
for a particular case* 

Bor the estimation of correlation eneigy in eqn. (5) 

“50 

we have used the expression of Nozieres and Pines-'^ in place 
of more lucent ones provided by Hedin et al. This is 
■!:ocause of better correspondence of with the observed 

value obtained by the use of the expression provided by the 
former expression. This aspect we will discuss in the next 
chapter along with the calculated results, 

Bor 'Calculating G-[n] for this case in eqn, (l) we have 
used the local density approximation and completely negle- 
cted the gradient teim. Incidentally, this is the general 
practice for evaluating the structure independent or one centre 
term in the pseudo potential perturbation theory where 
deviation from uniform density resulting in a periodic compo- 
nent is assumed to be small and so treated by linear response 
theory. This deficiency can be made up by treatrjag this 
deviation along with 'the uniform component together by the 
use of the density functional formalism given by the eqn,.(l) 
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where the gradient tern also oan be taken into account. 

One should also remember that this expression is formally 
exact, 3o the next obvious step is to try to incorporate 
such deviation in the parametric density functional forma- 
lism to evaluate both one centre and two centre terns accura- 
tely* 


(b) The periodic density variation : In place of a uniform 
density here we take a density distribution in pure metals 
incorporating a periodic component also. The form of n(^ 
should be such that it is continuous throughout the solid 
with no discontinuity in elope which will be unphysical in 
the sense that the higher order gradient terns will diverge. 
The form of n(r) which we have chosen for our calculation is 
given within a single WS all centred at r = 0 by 

n(r) = n^ + ^ cos k^r 

Since we req.uire that the function and its slopes are con- 
tinuous across the Wigner-Seitz cell boundary of r, we have 


dn 

dr 


^ 1 = 0 

r=H 


or 




lit 

1 


(4) 


when R is radius of the W.S cell. 
Thus, n(r) =11^+ 


( 5 ) 
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Tht constraint on n(r) is that the to'''al charge must 
be conseived i*e, the system is neutral electrically. So 
one has in each unit all 

n(r) ?r = Z (6) 

0 

Thus, if n periodic density component B^*s are assumed to 
be non-zero » there are (n+2) parameters 
which one is fixed by the constraint eq.n« (6). Since 
strong density inhomo genities are disfavoured (e.g. via the 
gradient tern), we assume for simplicity that only one , 
that corresponding to the smallest value 1 = I is non-zero. 

In this case the eqn. (6) gives 

_ ■ .. SP— ^ 2l£ll. cos IcE — 2 sin kBj (7) 

where B - and k = 

Now, from eqn. (4) we get kR = n and putting it in eqn. 

(7) n is evaluated leaving R and B as two free-parameters 
to be determined variationally . 

With the above two models for calculating the ground 
state energy of pure components we pass on to the problem 
of constructing models for alloys in the same spirit. 
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3*2 Step-Model 

In construoting this model we start with the spherical 
atoms of components with uniform density - "both solute and 
the Solvent, Prom the solvent matrix we remove an atom 
from the lattice site at origin and put the solute atom 
there. Pig, 3*l(3') shows the solute atom electiron density 
to be unifoimly n^ and its radius Pig. 3.1(b) shows 

the prepared solvent matrix with the atom at origin removed. 
The density here is also uniformly n^ and its atomic radius 
is Eg, Here the valency of the solute is taken to be lower 
than that of the matrix. Now if the solute atom election 
distribution is either expanded or compressed to fit into the 
rigid solvent matrix we arrive at a configuration of the alloy 
whose densi-ty distribution is shown in Pig. 3,l(c), The 
energy expended in this process is called the size factor 
energy contribution to the alloying. Still now no electro- 
chemical interaction between the solvent and the solute has 
taken place. Now if we start moving the step with two res- 
trictions namely (a) the electron density of the solvent 
matrix far away from the solute remains the same as that of 
the pure component and (b) the charge of the total system 
is conserved, we can arrive at a configuration given by Pig. 
3.1(d) which minimises the eneigy expressed by eq.n, (1), 

The energy involved in passing from configurations 
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(n) SOLUTE (b) SOLVENT 



(c) ALLOY (d) ALLOY 


Q -J,! t BOUBMATIO DUOiLfcM WB THB OHAROB 

(a) PURE SOLtJTB (A) IH 1 W3 OHM (b) ^R3 MLVBJT U>J[p* 

THH H3 OHLL AT/ORIOIK RBHOVBB («) 2?5?S5ii^®* 

PITS BIAOTLT US SOUTBST (d) AIlOT FOB OTBOT MIBIMull 

OOHPIOUR. '‘lOB, 
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correspoi-ding to Pig. 3.1 (o) to 3.1(d) is tlie electro- 
chenical oontribution to th.e alloying process resulting 
in a dipole at the surface between the solute atom and the 
matrix. This is accompanied by a lowering in eneigy when 
alloying is favorable. The density in this model can be 
expressed by 

n(J) = n^ ©(rg^ - r) n° (8) 


where is the distance from origin of the step and n^ can 
be found out from the charge conservation as 

.0 


= “b + 


3(2^ - 




(9) 


where S. and Z-r, are the valencies of the solute and solvent 

A Jd 

respectively. Now the densiiy is completely defined in 
terms of the parameter rg^. Following eqn. (49) of 
Chapter 2 we can proceed to evaluate different terms for 
finding out the heat of solution. 

The pseudo potential contiubution, TPl, .to 
embodied in the first teim of eqn, (49) of Chapter 2 is 
given in terms of Tg^ by 


^st - Eg 


3Z. rf 1 


r3 

^t 


R. 


] 

■B * 


( 10 ) 
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where A Z = and r^ is tlie core radius of Ashcroft 

local pseudopotential of A af fiie origin given by 
v^(r) .= 0 for r < r^^ 

= - ~ for r > r^ (U) 


TOl ='^. 


The electron-el eoti*on interaction term in the expression 
g^loy luy 2)01 is given by 

1^2 5 34 ^ ^st)-t-4'^n, 

^st + ^l-~b + T “ T" ^• 

- ^lt)J 

A Zid ^7 ^ ^ 

rr n 1^-2 C « rt . a . a Z, 

^ ~T- 


ji 


f 


10 


•St 


(12) 


+ 1.5 ZgizC — +^-^) 

® -5E" 5vl^ 




for r„ 


St ^ ®B 

The kinetic, exchange and correlation contribution to 
in the local density appiXDximation is given respectively by, 


IK = I5 4/3 + ^ 5/3 

IZ = -0.75(j)l/^ (Vg - 73 ^)] (13) 

“a ^o(n^) + (Vg - 73^) ^^(nj) 

47t r|^ ■ 47t R| 

where Vg^ = — j-2-=- and is gi'ven by the 

Pines and Uozieres expression for the correlation eneigy 


discussed in Chapter 2 
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Now taking the expressions given by eqns, (10), (12) 
and (13) together we get as a function of rg^ 


5^1oy = ipi + ici + TE + T2 + TO 
A 


(14) 


Now subtracting from it the ground state energy per atom of 
solute obtained from the uniform density grotmd state 
energy discussed in the earlier section of this chapter 
the heat of solution of the alloy per atom of A which 
can be minimised with respect to Here one should 

note that the matrix electron density is taken to be the 
same as that of the pure case, 

/■^ \ o . 

The assumption that as r -*■ » , n^r; ^ n^, i,e, 

the bulk density of the matrix B remains the same as in the 

pure B, needs to be examined. If it is not so and lim n(r) 

r 

= n^ then. A n = n^ - will be proportional to the concen- 
tration c of the solute in the alloy (from the charge neutra- 
lity). The atomic radius of the solvent in the alloy will 
also change (since each WS cell is neutral) leading to an 
increase in energy proportional to An^ i,e, o. But since 
the size of WS cell to which the solute has to fit in the 
alloy has changed, there will be an accompanying change in 
energy proportional to A n, i,e, c but this change wiJJL 
take place at o-oentres of solute making the change in this 
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2 

energy te-m proportional to c • Hence the total change in 
energy will be decided by a balance between these two teims 
and will be proportional to c and so in the dilute alloy 
limit one can see that since c 0 and so 0 

and the assumption that Lim n(?) =± nS is quite valid. 

Now one can argue that the neglect of 6n^ i.e, 
assuming uniform charge density of the matrix, may be quite 
serious because the effect of alloying will also result in 
modifying leading to an energy effect, fo estimate this 
energy change we examine the equation (l,3l) -vdiiGh shows 
that the energy effect will be of . second order in the 
chan ge of 6n^, Since the change in dn^ due to alloying will 
be quite small (excepting at the nearest neighbor sites 
which is taken care of by the construction of this 

energy change can be neglected. But if the concentration is 
high these arguments become invalid, 

3, 3 The Continuous Model 

The step model suffers from a serious handicap since it 
totally neglects the gradient term in the density functional 
eapansion. This term is by no means negligible and cannot be 
taken care of by a slight modification of density as has 
been done by Hodges and Stott^ (described in Chapter 1), On 
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the other haad the step model allows for a transparent 
physioal picture leading to separation of the size factor 
and the electro-cheiaical contribution in the alloys. 


To choose a parametric function for the purpose of mode- 

1 ling the charge density distribution in an alloy we looked 

20 

into the choice of Smith to model the charge density at the 
metal surfaces. Smiths form is exponential, i,e. it is given 
by 


n(2) 








for Z < 0 (in metal) 
for Z > 0 (outside) 


(15) 


where n is deteimined by charge conservation and p is the 
free parameter to be deteimined variationally , The resulting 
charge distribution for different metals is veiy similar to 
a tangent hyperbolic functiono The calculation of Lang 
and Eohn^ with a jellium model for the metal-surface also 
yields a charge distribution in agreement with the form 
obtained by Smith. In addition, the later calculation produ- 
ced small Priedel oscillations in charge density for Z < 0 
superimposed over the tang ert hyperbolic kind of distribution. 
In the studies of Bobinson and de-Ohatel^ also one finds 
a similar charge distribution for the jellium model Eohn- 
Sham oalcjulation in heterovalent alloys. From all these 
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C"J.C3ula'tioiis 11; is obvious tiiat a choice of "tangeii'b-’ 
byperholic functions is q,uii;e realistic for our purpose. 

So in tbis laodel the choice of 6n^(r) is 210 longer a 
step but a continuous tangent hyperbolic function and the 
electron density of the alloy systean can be written as 


n(r) = ng + 

B[l-tanh p(: 

r-ro)] 

+ I 6nS(^-Ei) 

(16) 




1?^ 0 


= HbCp 





= “b 

- n^ ©(Eg - 

1r|) 

+ 6n2(r-ST ) 

Wo ® ^ 

(17) 



l?l) 

+ B[l-tanh p(r-r^)3 

(18) 


B is a constant detemined by conservation of charge, p gives 
the elope and r^ the point of inflection of this distribu- 
tion. The last two parameters, p and r^, are the variables 
with resp3ct to which the grcjind state energy is minimised. 

Using the expression for given by eqn, (18) we 

can evaluate A H^from eqn, (2,49). Here, the two centre terms 
are no longer zero although in their evaluation also the con- 
tribution of 6n^(r) has been neglected because it is quite 
small. 

3«3.^ Details of Computation : 

For a given p and r^ we can find out B numerically from 
the charge conservation criterion yielding 
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B 


"a - 23 


U9) 


[l-tanB pCr-r^)] dr 
Tlae evaluation of B oociplotely defines n(r) and lienoe n-n (r) 

^ -LO O 

and n^Cr). The integration involved in eqn. (19) has to 
be done extremely acourately otherwise incomplete cancella- 
tion of long range teims in density functional foimalism 
will give rise to residual long range interaction giving rise 
to spurious results. In our calculation after deteimining B 
we have integrated n 2 ^^^(r) in the whole space and found it 
equal to with an accixracy of +.0001, To achieve this the 
above integral has been evaluated analytically from 3 ^ tq^and 
numerically from 0 to r^ by 15-point Gaussion quadrature 
formula. 


Now with the and n(r) obtained by the above pro- 

cedure we proceed to evaluate all the integrrls involved in 
eqn, (2,42) again by 15 -point Gaussian quadrature formula 
which we found accurate enough for our purpose by comparing 
it with the results of higher point quadrature formulas. 

The two centre integrals involved in eqn, (2,42) have been 
simplified analytically as given in CJhapter 2 leading to 
eqn. ( 2 , 49 ) and then these are converted into one centre 
integrals with the help of lowdin's a-function expansion. 
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The £H^(§,r^) thus evaliated is miulni^ed with the 
help of Powel's method of unco ns trained minimisation with 
two variables. 

To check the accuracy of the numerical procedures invol- 
ved we took a n(r) corresponding to the step model in this 
scheme and evaluated the energy and compared the results 
with the analytical results of the step model. The agreement 
is ezcellent upto the fourth decimal place giving us the 
necessaiy confidence in our computer program and the numeri- 
cal methods followed. 


Chapter 4 

Parametrio Variational Results and Discussion 

In this chapter ve present the results of our calculation 

with the paranetric variational scheme described in the last 
chapter. Both the pure metal and alloy ground state energy 
have been evaluated to determine the properties of physical 
interest. 

In Section 1 we present the general scheme of investiga- 
tion we have adopted, to achieve our primary aim to investi- 
gate the systematic trends in alloying behaviour. This will 
familiarise the reader with the fact that our purpose of 
the studies on pure metals emerges from the necessity of 
the studies on alloys. 

In Section 2 we present the cohesive energy, equili- 
brium atomic radius and the charge density for elements with 
different core radii in both the unifom gas model and 
the cosine model in the density functional formalism. We 
have also presented here our results on the estimate of the 
contribution to the cohesive energy of the periodic component 
of oharge density in the linear response theory to compare 
it with density functional results. 

In Section 3 we present the heat of solution of the 
alloys broadly under two categories (i) fixed matrix and 
the varying impurity and (ii) the fixed impurity and the' 
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varying natrix. Then an attempt has been made to compare 
the calonlated results with the experimental results. 

In this section we also examine the prevalent ideas 
on the effect of size factor and electrochemical factor 
in alloys. The discussion of the latter effect necessitated 
a brief discussion on chemical potential. 

In Section 4 we discuss our failures in the case of 
alloys with solutes of higher valency compared with the 
valency of the matrix. Our further attempts to form alloys 
with the solute having the gimund state charge density given 
by the cosine model are described. 

4.1 Scheme of this Investigation 

In the Ashcroft empty core model potential any ele- 
ment of a given valency can be characterised by its core 
radius, lilt hough this core radius is an artifact of the 
model having no rigorous physical standing, it is generally 
identified with the ionic radius of the elements. This 
can be understood in view of the cancellation theoran, i.e, 
within the span of this ionic radius the core state wave 
functions exist and, due to their orthogonality to the 
valence state, cancels, the effect of the ionic potential 
within this radius. If one calculates the r^, the radiv^ 
of the sphere containing one electron in the element, from 



11 ^ 

■bhe rainiELG-l property of the ground state eneigy by a uniform 
electron gas model one gets an approzimately linear curve* 
This linearity remains intact even after assuming a linear 
response screening excepting for very low value of the core 
radius. So ^ in general one can take a systematic variation 
of r^ to be equivalent to a regular variation of r^. Now the 
atomic radius being directly proportional to the r„ for a 

D 

given valency, r^ becomes a measure of the atomic size of the 

elements. For the same valency r is a measure of the 

o 

strength of the potential also* Hence a systematic variation 

a 

of r^ is equivalent to ^variation of the strength of the ionic 
potentials as well as the charge density for the same system 
of electron gas as we go along a homovalent series, Inciden 
tally, Inglesfield has taken the difference in the depth of 
the pseudopotentials in homovalent binaiy alloys as a 
measure of difference in electronegativity. So by vaiying 
this parameter r we are able to systematically study the 
effect of size-difference and electro chemical factors in 
alloys as it will be detailed in later sections, 

- Bqn, (3.2) of the last section expresses the ground 

state energy B as a function of R which in turn is pro- 
g ^ 

portional to r_ for a given valency. This B will also be 
a function of r which is characteristic of an element 
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and hence xor a given we Ccja nininise B with respect to 
Tg to get the Tg vg rs curve for a given valency. Singh and 
Young have already calculated such a curve including the 
contribution of the periodic conponents of density in linear 
response theory. ¥e also calculate the sane results essen- 
tially to check the accuracy of our nunerical procedure. But 
we use No zieres— Pines expression for the correlation esneigy 
and Hubbard— Shan dielectric function in place of the Hartree 
dielectric function. Our results agree in general with that 
of Singh Young including the non-linear nature of this 
r^-rg curve at low r^ as shown in Pig, 4.1, Since our nodel 
for the alloy conpletely neglects the periodic oonponent we 
nininise the eq.n. (3*2) i.e, the energy of a sphere of 
electron gas of Tjuiifom density and obtain the r^ vs Pg curve 
as shown in Fig. 4,2. Ihis curve is almost linear providing 
a good basis for the study of alloying behaviour along these 
curves. 

The non-linearity in Fig. 4.1 probably indicates that for 

the _ 

such low value of r^ pseudopotential becomes strong enough 
to nake its treatment by second order perturbation theory 
invalid. Otherwise the general trend of both Figs, 4.1 
and 4.2 are q.uite understandable in the sense that for the 
sane r^ the higher the valency | the lower is the Tgj 
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beoause hi,‘?;her valency gives stronger potential to pnll the 

electrons towards the ions and hence results in higher den— 

sity of electrons i*e« lower r « 

s 

One should note that these results given in Pig. 4*2 
differ fron that of Ashcroft and Lang retl^ Although the 
agreement of cohesive energy for a given r„ and valency is 

o 

better here as it will be discussed in the next section. 

4.2 Cohesive Energy of Pure Metals 

As described in Chapter the cohesive energy in pure 
metals has been studied in two different models (a) uniform 
electron gas model (b) the cosine density model. The moti- 
vation for talcing resort to the latter model is also descri- 
bed. Here we will present the results obtained in these two 
models in the following two sub-sections, 

(a) Uniform electron gas model : 

To give an idea of the nature of variation of the 
energy in this model with respect to r^ we p 2 ?esent in Pig* 

4.3 S' sample of our calculation for Z = 2.0 and r^ = 1,0, 1,7 
respectively* Also we show the contribution of the periodic 
variation of density in the linear response theory. The 
curves denoted by U represent the structure independent 
energy or the energy in the pure unifoim electron gas model* 



Cohesive energy in a.u. 




rs In Q.u. 


PliJ 4*3 * VAillATlON 0^ MBHQT WITH Tbi O-UKIFOflM JSIiBOTHOH QiS 
KODliLj K-I*^SJUa HBSPOKSa THBOiff POH DIPFJSRBHT 
VALBHgiaS (a) t • 1,0 (b) 2 * 2*0 ISO (o) Z » 3,0, 



l‘6 ' 2.0 3-0 4.0 


rs In a.u. 

no 4,4 I TTAHIATIOH OP OOHBOITB KNfiHOT WITH re* (D U-.lTi'{IP01<^: 

. BLBOTHON 0X3 MOPfiL# (li) LINBAft R.:jilF0i-':5ii THPa:; 

POa DIFPE'HKNT VALSNaiJSS WITH OOHRiUTiON 

MBKOT 0IVjSi\ BI 0MXi»'30ir3X XH (a) Z » 1.0, (B) ^ 2<-0| 

and (o) 3 - 3.0, 
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Fig. 4.4 s'lows the variation of cohesi-ve energy i.e, the 

position of the nininun of Fig, 4.3, with respect to r 

s 

for the unifom gas model by the firm lines for different 
valencies. The observed value of r^ and cohesive energy of 
Some of the simple metals are also shown in this figure by 
circled dots for a comparison with the calculated curve. 

For Z = 1,0 we could not show the contribution of the perio- 
dic variation of density because it is so small that it can 
not be shown separately in this scale of the figure 4.4. 

There are several expressions available for the corre- 
lation energy of the unifom electron gas. We have used the 
expression given by ITozieroesnd Pines (in atomic units) as 

— “0*^575 + 0.0155 In Xg 

per « electron. We also tried to use the interpolation 
^ 29 

formula used by Ounnarson et al and the results for Z = 2,0 
are plotted in this Fig. 4.4 as the dash-dotted curve marked 
by L, This curve shows similar results as the NP exchange but 
the agreement with the observed cohesive energy is inferior 
in this case* So we abandon this interpolation fomula and use 
the one given by Uozieres and Pines consistently throughout 
this study . 
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Tiie r,greenent of the observed values of cohesive eueigy 
and r_ with the calculated curve is q,uite excellent conpared 

P 

to the previous works of Ashcroft'^ et al and that of 
D, Wearie^^. The foraer study did not achieve such good 
agreenent in cohesive energy whereas the latter study starting 
with the Aninal'U— Heine Model potential for differoat elenents, 
did not get good agreenent with the observed values, 

Por the purpose of our study if we define the core- 

radius r of an elenent as the one which corresponds to the • 
c 

observed atonic radius in the unifom electron gas nodel, we 
find a disagreenent with the core radii derived fron the ob- 
served resistivity or the femi-surface data. This is expe- 
cted on two counts: (i) due to ignoring the periodic contri- 
bution which has been taken care of within the franework 
of the second order perturbation theoiy while deriving the 
core radii fron the FS or resistivity data (ii) due to the 
fact that FS or resistivity data reflects the effect of the 
ionic potential on the states near the FS and so will sim- 
late the scattering of the electrons in these states only, 
whereas the core radii found out on the basis of cohesive 
energy will give a better nodel potential when the scatter- 
ing of all the electronic states in a band fron the ionic 
potential is inportant. 
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If we define r as indicated above then our calCTiLation 

yields a nuch. better agreenent with, the observed value of 

E (corresponding to an atonic radius) than any other 
S 

previous study. 

(b) The Periodic Density Variation : 

The unifoin charge density nodeli although is well known 
to account for nost of the cohesive energy ( 90-95 percent 

by earlier accounts, our result shows it to be even higher) 
but still a snail contribution fron the periodic oonponent 
of density is observed# Since we are interested finally 
in the heat of solution of the alloys which are q.uite 
snail (in cases nearly of the order of the cohesive energy 
contribution of the periodic conponent) it is desirable at 
least to understand this contribution if not acconodate it 
in our study of the alloy. Oustonarily , this has been 
studied under the linear response theoiy but we find it 
quite unreliable in the linits when of snail r^ ahd large r^^ 
coupled with the enpty core nodel pseudopot ential. So, 
we study the ground state energy with the fomalisn of the 
density variational nethod along with, the gradient tern a 
described in Chapter 2, 

(i) The results of the linear response theory J 

The curves denoted by B in Fig. 4.3 show the variatrcn 



82 


of Energy ■with r^ including the contribution of the periodic 

variation of density for Z = 2,0 and the core radii r_ = 1.0 

and 1.7 for fee structure. The calculation for the h#p 

structure shows alnost the saoe values of E for a given r^. 

It is expected because the energy involved in alio tropic 

change between these two structures are observed to be eztrene— 

ly saall if the volune per aton remains the sane. The 

behaviour of these B vs rg curves shows that as r^ becomes 

larger the total energy for a given rg is higher because the 

attractive ionic potential becomes weak. Also for a given 

Z as r increases the energy contribution of the periodic 
c 

component of the charge density increases. When the core 
radius is large, only by piling up electrons in the attrative 
region outside the core can the systems reduce its energy 
giving rise to a g 2 ?ound state density with strong spatial 
dependence i.e, periodic component is large. This is obser- 
ved as a systematic trend in Fig. 4.4. Here the dotted lines 
give: the energy including the periodic density contribution 
at the r„ obtained by minimising the energy of the uniform 
density electron gas. Iliia procedure is justifiable on the 
ground that the periodic density contribution is^sxtrenely 
snail compared to that of the unifom conponent^or high r^ 
the periodic oohtributioh increases significantly and the 
above procedure nay not hold. This increase iu the 
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pBriodic contribu’tioii with -loy > howBver, ha an artifa,ct 
of the ffi'-ipty core pseudo potential. 

This periodic contribution increases with valency also, 
as expected, because the periodic potential gets stronger 
with Z. But when Z is large the potential nay be strong 
enough to the extent that the linear response theory if 
applied nay over estimate the periodic contribution. 

In general, it appears that the linear response theoiy 
coupled with the Ashoroft empty core pseudo potential tends to 
exaggerate the periodic contribution increasingly 
with increasing r^ and for the low r^ and high valency it 
gives a wrong estimate (Fig, 4.1). 

From Fig, 4.3 one can also see that the rate of variation 
of this periodic contribution with r„ can dominate over the 
rate of variation of the structure independent unifoirm ele- 
ctron gas energy with r_, for higher r. 's and as a result 
one can find an r_ corresponding to the minimum in total 
energy quite different from that corresponding to the minimum 
of the uniform electron gas eneigy. Since for the reasons 
mentioned in the above paragraph^ we suspect the estimate of 
the periodic contribution under the model of empty core poten- 
tial, we took the model of \anifora electron gas without the 
periodic variation of density to depict the electron 
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distribution of pure metals in the context of alloys. Foi^ 
tunately also it gives quite a good agreement vd.th the 
observed cohesive energies justifying to an extent this 
model. This model does iwxLty in the cases of alloys 

with impurities of higher valence which will be discussed 
in the appropriate place. As a result it becomes important 
to deviate from this model for pure metals in an effort to 
find a satisfactory treatment of the above alloys. In the 
next section we incorporate a periodic variation in densi-ty 
and evaluate the cohesive energy of pure metals by the 
density fimotional formalism in place of linear response 
theory. 

(ii) The Density Punctional EesiiLts ; 

In this scheme we use a density distribution of the 

form n + B cos kr with three pairameters n , B and k as 
0 

described in the last chapter, n^ is determined by charge 
conservation and the other two parameters are determined 
variationally usirg the minimal property of the ground state 
energy. The numerical scheme is similar to that of "the 
tangent hyperbolic density model of the alloy. In this 
scheme we find out the dominating contribution from the set 
of k’s which obviously will be the one having half the 
wavelength within the atomic radius. Other smaller wavelengths 
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will have larger gradient tena contribution reducing their 

effect on cohesive enei^y. This dominating k is related to 
atonic radius R . by 

-t 

R - it/R R, (= 71 /Ic) and B are detemined by nininun 

B 

Fig. 4.5 shows a plot of R, . as detemined by this method ^ 

for valencies Z = l.o, 2,0 and 5,0 respectively. From this 

figure we can see that at core radius (r^) greater than 2.1 

this R^^ Vg r^ curves deviate' from the linearity observed 

for lower values of r^. The dotted curve at the top shows the 

^c 0^ the uniform electron gas model for valency 

Z = 5 . 0 , The comparison of the dotted curve with the 

Z = 5.0 curve for the model shows that the of R 

increases with increasing r ( 1,5 percent at r =0,8 and 17 / 

percent at rc=2.7 - ® 

This may be explained in terms of the periodic effect we see 

for the empty core potential. But this deviation is quite 

small with low r^ and increases with the increase in valency. 

Figure 4.6 shows the difference A E in cohesive energy 

between the uniform gas model and the cosine model as a 

function of r^. We see that there is a general increase 

in the periodic contribution with increasing r • TiiTf:! 

o 

can be understood in view of our argument that if r is 
large electrons have to pile up in the small attractive 



Atomic radius(R) in q.u. 



P 10 4.5 I VAaiATION OP BQUILLIBaiaM ATOMIO HADIUS (R) VlIH 

OOHS RADIUS (tq) in Ea 003 INB MODRL FOR VAIJWOIja 

(a) Z » 1.0, (b) Z - 2.0 AND (o) Z - 3.0, (d) 

FOR UNE'ORM BLBOTRON OAS MODUD 1 -« 3.0. 



P 10 4.6 t VARIATION OF DiFFiiRENOB IN OBOUND fci'JJATB aWSROIBS 
(AB) NITE OOitS M0IU3 (rg) BSTVBISN TilB 008110! 
MODBL AND THB UtllFOM ELBOTHDN OAS KDDBl FOR 
(a) Z » 1.0, (b) I « 2.0 AND (a) % m J.O. 
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region r^ < r^ < B- » end hence "tiie periodic component and 
its contribution to the cohesive eneigy increase Pig. 4»7 
shows the variation of AB with the valency Z, AE increases 
almost as for r^ = 1,7. If we examine in the same 
figure the variation of -nQ/B with Z we find that it satura- 
tes at high Z( > 3»0)f i,e, with Z = 3*0 and b^ond B 
increases in proportion to n^, Big, 4»8 shows the variation 
of -%/1b I with r^ for different 1 , starting from r^ = 1,2. 

We find that with increasing r^ the periodic component 
(B| increases faster than n^ (i.e.lBj/n^ increases until it 
reaches a certain value then again lB(/n^ decreases 
(probably dis continuously ) . This is not expected, since the 
g,s,e. E. is a smooth fimction of r ,B, and B, Most likely 
this effect is due to an incomplete search for the true 
minimum, i.e, the minimum obtained may be a relative minimum 
but not the absolute or global 'Minimum. We have not 
pursued this further since the' energy differences are small, 
and the values obtained here (beyond the jump) are not used 
in further computations. For small r^ = 0.8, 0.9 we find 
that B is positive indicating a pile up of cjhaige at the 
centre. This is because the strongly attractive region of 
the potential is too near the core leading to a pile up of 
charge there and this continues even inside the core because 
gradient term does not allow the charge density to fall 



as 
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sharply. The cx)Sine fora of charge density nay not have 
correctly represented the density in this case. 

If one conpares the periodic contribution to the cohesive 

energy fron the oosine nodel and that fron the linear response 

theory one can easily see the tendency of linear response 

theory to overestinate the periodic contribution to B in 

S 

general, Bor low rQ(= 0,8 or 0.9) we find (n^/B) positive 

and large indicating a snail pile up of charge at the core 

region. The linear response results in this region shows 

a deviation fron its nomal behaviour as shown in Big, 4.1 

and since the potential is quite strong for low r these 

c 

resxilts are suspect, 

4.3 -The Heat of Solution of Alloys in the Paranetric 
Variational Model 

The calculations have been done on the following systms: 
fized natrix with varying alloy elenent, and fixed alloy 
elenent with varying natrix. By fixed, we nean that the 
valency 2 is fixed as well as the core radius r^ (and thus 
the atonic size). By varying, we nean that for a given Z, 
the core radius r^ is varied. Thus one is able to study for 
a given valence difference, the effect of size difference 
(e.g, difference between the equilibyiun atonic radii ^ 
and Bg) as well as the effect of the chenical potential or 
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elec’bro—EL^aga'bi'vi'ty difference. The lafter depends on "tiie 

electron density wliiGh. varies snoothly with r. Such a set 

o 

of calculations should bring out clearly any systenatic 
trends in the alloying energies etc. The observed heats of 
solution for different alloys correspond to the heats of 
solution for particular values of r ’s for a given Z. - Z-p 
conbination. Even though the actual values of the heht of 
solution nay not agree well with experinent, the trends 
observed nay have deeper validity to be confimed by ezperi— 
nents, synbolically , the natriz-impurity conbination can be 
represented as MZ^ - foiner set (fixed natrix varying 

inpurity) contains results for the following: = 5 and = 1,2 

Z ^ 2 and = 1,3 and Zg = 1 and = 2 and 3* In the latter 
set (fixed inpurity vaiying natrix) we have the results for 
Zg = 2 and = 1, 

The results can clearly be classified according to the 
relative valencies of the natrix and the inpurity, i*ei 
Zb > Zb = and Zb < Z^., The variational paraneters we 
have used to describe n(r) are such that only in the last 
case can we study the alloying effect. The reason is as 
follows. In the case Zb > Z^, we generally find no charge 
transfer, but only a redistribution of n(r) within the WS 
cell corresponding to the impurity. Ihat is, in the step 
model the step occurs at r^^ < ilgr and in the tank model 
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for n(r) tlie turnover point r^ < R^, This nciification of 
the ground state profile might be due to the host environ- 
ment, or it might be intrinsic, i»e, present in the pure 
metal A aS shown in Section 4. 2, Because of the lack of charge 
transfer, this does not fall into the canonical alloying 
scheme, we however, discuss the results for this case in a 
latt.or section in some detail. In the c^se = Z^, it can 
be shown, that one can have no step in n(r) nor can one have 
an additional tanh kind of distribution. This follows 
simply from .charge neutrality. In our variational parametri 
zation* scheme the only contribution to is due to di- 
fference in atonic volumes, Olearly, a more sophisticated 
density profile is needed for the homovalent case, 

¥e thus turn to the case Zg < Z^, where we find charge 

transfer and thus alloying in the conventional sense. We 
present our result as follows i 

(a) Fixed matrix with varying alloy elements ; 

Figs. 4.9(a), (b) and (c) present.- respectively the 
results for the matrices with valencies Zg = 1,0, 2,0 
and 3.0, the core radii in these three cases have been 

r® 1.12. 1.2 and 1.2. Apparontly.ihaBG . curves represent 

c ' 

widely varying trends for different - I2»^ oonbination. 
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But in most of the casea the results of the atep model 
follows the genenal trend of the restilts of the continuous 
tangent hyperbolic model although quantitatively the 
results of these two models are different* Olearly it 
indicates the importance of the contribution of the gra- 
dient terms in alloying. The values of r and r are also 
differsnij because tlie "bauh model does not merely relax tliQ 
step keeping its turnover point r^ at r = r^^ but it shifts 
the turnover point also from r = rg^ to r^ < rg^ in all the 
cases (See Table 4*1)* In the M5-I1 and M3-I2 set the heat 
of solution of the step model is greater than that evaluated 
by the continuous model but for M2-I1 we get the contrary 
result' • As we go from the former model to the latter we 
find that two opposing tendencies develop: (i) changing the 
charge transfer contribution to the energy, i,e* the 
negative energy contribution of the surface dipole energy 
and (ii) the positive energy contribution from the gradient 
terms. As a result the charge distribution n(r) readjusts 
itself to arrive at a new minimum* 

The widely varying trends of the heat of solution 
aSj. Vs the core radius of the solute r^ as given in these 
figures are due to the fact that for the sets where 

< 2* thet^ are no eleotro-ohemioal interaction between 
the solute atom and the matrix beoeuse it was found that 
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and so also r^. This is equivalent to a redis- 
tribution of charge within the impurity cell which since 
deviates from its ground state configuration should give an 
increase in energy* But we find that in such cases also 
there is a negative region of AHj which is puzzling* Our 
attempt to investigate this problem will be reported at the 
end of this chapter* If we ignore these results for Zg < 
the rest of the sets (i*e* the results for > Z^) follows 
a definite trend which we plot separately in the nezt 
figure* 

In Fig* 4*10 we show a plot of the atomic radius 
of the alloying elements VS the heat of solution AH^ at 
infinite dilution for these cases where alloying has taken 
place. It includes three sets of variational results 
denoted by and M2 - H- This figure also 

includes one non-variational result for M3 - 13 for the 
pu2?pose of comparison. 

These results show that as an alloying element with 
differing lower valency is introduced in a matrix the 
AH-p VSR, curve becomes aS 3 niimetric. This clearly is due 
to the electro -chemical effect because in the alloying with 
components of same valency the size factor only has 

and leads to a nearly symmetric AHj VS 


got a role to play 




?xa 4aO * ?AHIATIOIi OF HBAT 0? SOiUIIOH, AIt. WIIH SOLUTB 
Aimia HAiJlUS, iQB. (a) M3 - 13, (b) M3 - 12, <o) M3 - U 
(d) M-i - U ST3T£M3. SXPKIlIMmAL RilSlJiTS* 
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curve as given by the M 3 - I3 set of Rig. 4.8. But whea 
we go to M 3 - 12 or M 3 - n we find curves for iC\ H ^ whose 
assymetry increases with increasing difference of the valen- 
cies /A Z of the constituents. The minimum of these curves 
shifts to a lower and lower size away from the atomic size of 
the matrix with increase in 4 \ 2 * In other wo2?ds the electro- 
chemical factor in these cases makes alloying with impurities 
of the smaller atomic radius ratio (radius ratio l) more 
favourable compared to that of the elements with larger atomic 
radius ratio (radius i*atio ^ 1) so, the nature of 15 percent 
rule of Hume— Rothery should incorporate an asymmetr7 due to 
electro-chemical factor of the non-homovalent alloys and this 
limit of 15 percent should increase with the difference in the 
valencies of the constituents. One notices in Rig. 4.9(c) a 
drooping of the M 3 — 12 and M 3 — H curves corresponding to 
the whole charge of lower valence impurity getting stuck to 
the charge of the matrix leaving an almost charge free region 
in the centre. This is equivalent to readjustment of charge 

1 

in the impurity cell and so should lead to a positive energy 
change contrary 'to that shown in Rig, 4.9(c) for M 3 ~ 12 for 
the step model but of course the continuous model rectifies 
this situation. This particular trend may be an artifact 
of this potential due to similar reasons as mentioned for the 


pure cas e • 
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4«H sliows Ijlid offset of olisjogo of (ooro 
radius of inalirii) on '"the^^j VS cuivo* A reduoiiion ^ y* 
results in a reduction of atomic sizie of the matrix and 
hence this curve shifts ,as expected to the lower size. 
There is no significant change in the depth of the minimvim 
which appears to be a strong function of the difference 
in valencsy of the constituents only, 

(b) ^ixed impurity with vaiying matrix ; 

Figure 4#12 shows theZSHj VS Eg results for a system 
with the matrix and the impurity having valencies 
2»g = 2,0 and 2^ = 1,0 respectively, the core- 

A 

radius of the impurity being r_ = 1,5, The corresponding 

V 

E^ is shown by the dotted line in the figure. This curve 
does not show any new aspect but only confirms the nature 
of the asymmetry observed earlier. The minimum always lies 
where ^ Eg and the alloying with solvents of larger 
atomic radii is more favourable th an that with the solvents 
of smaller atomic radii. 

The maximimi calculated negative value of the heat 
of Solution is in the range -0,02 to -0,03 a,u, or 12 to 18 
ECal/mole, This is also in the range of maximum observed 
value. This number is very small in' comparison to the 
cohesive energy ( 5 percent). Often, observed heats of .* 
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Solution are smaller* Thus one can expect that there will 
be many oases of disagreement between theoiy and experiment* 

Since we have not calculated the heat of solution 
except in g for the Al— base and Zn— base alloys we devise a 
method for extrapolating the heat of solution for other 
systems* Since we see in Pig, 4*11 that the effect of chang- 
ing r^ of the matrix produces a shift of the Y3 curve 
keeping their shape almost uinaltered, we plot all the data 
of A Hj against Rg - in Pigs, 4.13 and 4*14^ In Pig* 4*13 
the squared points are the data for M2 - H systems with 
varying matrix and fixed impurity* The triangular and cir- 
cled points are from the fixed matrix and varying impurity 
set for the M2 - H systems with r® = 1,0643 and 1,2 respe- 
ctively, When the size mismatch is small all the points 
fall on a single line. But they start diverging when the 
contribution of the size factor to the heat of solution 
becomes significant* Pig* 4*14 presents similar curves of 
A 73 ( % - ) ^or the M3 - 12 and M3 - D- systems, 

( c)0omp arisen with Experiment • 
i) Homovalent Oase i 

Before discussing alloying in the heterovalent oase, 
we collect here e:^erimental results for homovalent sinqple 
metallic alloys. Here, within the parameter space 
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chosen 't 7 us for the variational method, there is no- 
step or tanh density profile. The only contribution to the 
heat of solution is from the size or atomic vol-ume difference, 
this being necessarily a positive teim. 

In Table 4.2 we show for various valencies, the heats 

.A 

of Solution for AB alloys (both dilute^^and dilute B limits), 
and the conditions under which they are measured. The sign 
of heat of solution is indicated if that is the only tiling 
known. Ve see that while in a number of cases, the heat 
of solution is positive, there are nearly as many examples of 
negative heats of soluliion. Thus, the size difference 
cannot be the only or dominant factor influencing AH Ve 
discuss below the various alloys. 

In the noble metal alloys, obviously the empty core 
description is inadequate. The. noblo nettilo have a just 
filled d shell, considerable s-d mixing, and the ionic 
radius is comparable to the atomic volume radius. Therefore, 
there may be electron overlap effects not described in a 
simple pseudo potential model. 

It is kno-wn that in Mg, because there are very- 

few core orbitals, the core orthogonaliaation is not very- 

effective and the ’pseudopotential* is q-uite strong, In- 

. » 

deed, many properties of Mg are not well described by a 
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TABIB 4,2 

Partial Molal Heat of Solution of Homovalent Alloys at infinite 
EXPERIMENTAL RESULTS dilution 


System In E cal/gm.atom 


Type 

A 


B 

5a 

5b 

TEMP, STATE 

Ml-n 

Ag 


Au 

-4.05 

-4.85 

800®E, Solid 


Ag 

- 

Ou 

9.00 

6.00 

1052°E,3olid 


Au 

- 

Ou 

-4.315 

-2.84 

720°E, Solid 


Na 

- 

E 

-we small 

CD 

CO 

liqTOid 

c\J 

M 

1 

CM 

Od 


Mg 

-4.13 

-2.2 

543°H,-' Solid 


Od 

- 

Hg 

-1.88 

-ve 

600°E, Liquid 


Od 

«- 

Zn 

-2.14 

2.25 

800°E, Liquid 


Hg 

-- 

Zn 

-we 

0.71 

608°E, Liquid 

M3-I3 

In 


T1 

.38 

.69. 

623°E, Liquid 
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pseudopotential. In case of Od and Eg, both group II B 
elements, factors si m ilar to those operating for noble 
metals (large core size, some s—d mixing) are present and 
thus can lead to negatived 3^, e.g. for OdHg, 

The systems E— Na, Od— Zn, Hg— Zn> InTl, in which one or 
both the constituents' are simple metals well described by 

A 

a simple pseudopotential, all have opositive heats of 
solution, ■ In the Na-E, Od-Zn and In-Tl systems the calcu- 
lated values ofAEjj^ = 1,25 = 5*56 and^ = ,707 

respectively in the unit of E.cal/gm, atom* These values 
are of the same order as the experimentally observed ones. 
The observed trend of increase between these values i,e. 

^ maintained, 

ii) Heterovalent Case (MZg - IZ^ with Zg')^ Z^) 

ii-a) M? - n : 

As an example of the series of systems studied by us, 

valence 1 and has a 

we discuss here the case ^ere the impurity or solute has a^t 

given atomic radius (and thus pseudopotential cjore 

radius r"^). The host has valence 2, and its equilibrium 
o 

atomic radius is varied by vaiying r®. This enables us to 
generate a class of M2 — H alloy systems with different 
size and electro— chemical potential differences. 
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In Table 4*3 we show experimental results for M2 - H 
systems* The sign is shown if it is the only inf o nnation 
available* The calculated results are also shown in the 
' table. Actually# the calculations have been carried out for 
2n based alloys. The numbers for other cases (Od and Eg 
based alloys) have been obtained on bodily shifting 'Hie Zn 
curve by an amount equal to the difference in atomic radii. 

This procedure gives reliable numbers, piovided the compressi- 
bilities or size factor effects of the S = 2 hosts do not 
differ widely and provided these do not dominate the heat of 
solution. We have shown earlier (Fig. 4.13 f Section 4. •3(3') 
above) that this is indeed the case for the typical example 
described there. 

We see that except for the case of Hg-K the sign of the 
heat of solution is correctly given. The right order of 
magnitude is found for the only two alloys for which Hj 
is known experimentally# i.e,0dTCu and Hg N’a. As will t)e 
seen from examples discussed below, we tend to get positive 
heats of solution for oases of Solute oompression, i.e. 

\ Experimentally the heats of solution in such 

cases tend to be more negative than the theoretical value. 

ii-b) M3 - 12 and M3 - II Systans s 

Table 4*4 shows the observed and calculated values 
for the above class of sy-StoaB. Actually, the 
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TABLE 4.5 

A Hj for M2 — H Systans 


1 

M2 

System 

n 

Bipt. 
lues AHj 
iu a.u. 
atom of 
state 

Oalculated 
values /jHj 
in au./atom 
of state 

in 

Temp. State 
for the ezpt. 
results 

Od 

Ag 

-ve 

-.016 

+.256 

67%Solid 


Am 

-ye 

-.016 

+.257 

700°E:, Solid 


Cu 

-,0025 

-.01 

+.597 

875°E,liquid 


Ea 

-ve 

-.008 

+.673 

Solid, 298°K: 

Hg 

Ag 

-ve 

-.014 

+.35 

Solid,292‘^E 


Au 

-ve 

-.014 

+.351 

Solid, 550®E 


K 

-.0398 

+.045 

-1.51 

liquid, 600°E 


Li 

-ve 

-.016 

+.094 

Solid, 292®E 


Na 

-.032 

-♦016 

-.ar9 

Liquid , 298.15^E 


Ag 

-ve 

-.016 

-.104 

Solid, 87 3^E 


Ax 

-ve 

-.016 

-.103 

Solid, 36 3°K: 
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TABLE 4*4 



IX 


for M3 - n 

and M3 - 

12 Systans 


Type 

System 

A.N HT(Bxpt) 4iSj(Qalcu- 


Temp, State for 

M 

T 

in a.u/ 

. lated 
in an,/ 

the ezpt re- 
sults 




X 

solute . 

solute 





atom 

atom 



M5-I1 

Al 

Ag 

-ve 

.032 

-1.018 

Solid, 820^E 



Cu 

-.024 

-.022 

C\J 

# 

-f 

Liq,1423°H 


In 

Ag 

-.0019 

-.02 

-K).47 

Liq,723^E 



An 

-.022 

-.02 

+.471 

Liq,723°K 



Ou 

-.0013 

-.0015 

+.811 

Liq,723°K 


T1 

Ag 

+.0055 

-.016 

+.575 

Liq,723°K 



Au 

-.00008 

-.016 

+.576 

Liq,723°E 



Na 

-.03117 

-.021 

-.354 

Liq,648°K 

M3- 12 

A1 

Zn 

+.0050 

-.014 

+.086 

Solid, 625^K 



Mg 

-.021 

+.01 

-.355 

Liq,723°H 


Ga 

Od 

+.0049 

-.008 

-.104 

Liq,700°K 


Zn 

+.0019 

-.011 

+.256 

Liq,723^E 


In 

Hg 

-.0028 

-.014 

+.120 

Liq,433°E 



Od 

-.0018 

-.013 

+. 214 

Liq,723 K 



Zn 

+.0042 

+.013 

+.574 

Liq,700 E 


T1 

Hg 

-ve 

-.012 

+.225 

Liq,293°K 



Od 

+.0028 

-.007 

+.319 

Liq,673 E 



Zn 

0067 

- 

+.679 

Liq,1099 E 
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oalcjulations were done for iU. based alloys -bhe numbers 

for Inj Ga and Tl based alloys have been obtained by shift— 
in^ the^ H j curves as described earlier above (This section 
4*2(ii))* Again, the sign of^Sj is generally given 
correctly* The heat of solution is often overestimated in 
many good pseudopotential oases where there sire no problems 

due to large core, s-d mixing, etc* the agreement is fairly 
good* In many cases, we obtain the correct trend, e*g* 
in In-based alloys, as the solute changes (Hg, Od, Zn), the 
size difference increases, and so does the heat of solution. 
Our obtaining too large negative values may be partly due to 
inaccurate estimate of the size factor effect. The trend 
in Tl(0d, Zn) is also correctly reproduced, 

(d) Factors Affecting the Alloys : 

In the first chapter we have already discussed the 
state of present understanding of the factors affecting 
the alloying process namely, the size factor effect and 
the electro— chemical effect* In Chapter 3 we have presen- 
ted an account of the variational models - Ihe step model 
and the continuous model* It has also been pointed ottt 
there how in the step model these two factors can be 
separately estimated and assessed inspite of the crudity 
of this model* In the continuous model this advantage 
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is lost^ In the following sub-sections we are going to 
discuss these two factors in view of our calculated results^ 

(i) Size Pactor Effect • It is, as ezplained in 
Chapter 3 through figure 3%l(o), the energy spent to p 2 ?e* 
pare a solute atom so that it fits snugly in the solvent 
lattice site* With the help of eq.n» (3«14) we can write 
this energy A Eg ^ as 

( 1 ) ' 

< ^St = % 

In our computation we tajce a monovalent solute of 
core radius r^ = 1.2 and put it in a matrix of valence 2 
or 3 (M2 or M3). Ib.e matrix is immaterial here since what 
we study essentially is the expansion and compression be- 
haviour of an atomic cell of the solute of valency 1 and 
core radius 1,2. The results are presented in the form of 
plot of ABgp VS (R^ - Rg)^ in Pig, 4.15. The points 
marked by 1 represent the path of energy during compre- 
ssion, ¥e also repeat the calculation with solute core 

radius r”^ = 1,5 to study its effect and these results are 
c 

identified by squared points in the same figure. 

The general features of these curves are as follows • 
(1) There is a linear ;region when the deformation is 
small (2) This linear region extends upto larger 
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2 

expansion as compared, "to ' < compression* 

( 3 ) liie compression involves higher energy than expansion 
for the same amount of defoimation, (4) As the core radius 

increases the strength of the material reduces and hence 

A A 

the r - 1.5 curve lies lower compared to that of r'r = 1*2, 

O o 

These characteristics described above are quite 
expected and can be understood easily. The point we shoiild 
note is the severe limitation of the elastic continuum 
theory to describe solid solution when Eg 

Eg this theory should work well to quite a large 
extent of size-^aismatch. The general belief of symmetric 
behaviour with respect to (H^ - Eg)^ is not confirmed. 

This explains the asymmetry in the M3 - 13 curve of Eig. 
4.10* This asymmetry is much more pronounced in 
curves of Pig. 4^10 for no n--homo valent system because of 
the additional effect of the electro-chemical factor, 

(ii) The Electro-chemical Pactor » The electro— chemical 
contribution to the heat of solution is the energy change 
connected with the charge transfer between the solute cell 
and the surrounting solvent matrix. In the step model when 
rgt = Eg i,e. the configuration where the solute atom has 
undergone a dilation or compression to fit the solvent site, 
there is no charge transfer and the energy to be expended 
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to attain this configui^tion is termed the si25e factor con- 
tribution as explained in the previous sub-section. It 
corresponds to Fig, 3.1(c). Now this charge density is not 
the groimd state charge density in the alloy and so a further 
re-distribution of charge takes place leading to a charge 
transfer. This results in a surface dipole at the interface 
between the solute and the solvent. This reduces the energy 
of the system. The charge transfer i,e, the dipole forma- 
tion, will be opposed by an increase in the kinetic and the 
electron-electron interaction energy of the system, A b a l an ce 
between these two opposing trends will lead to a minimum 
corresponding to the groTond state energy of the alloy (Fig. 
3,l(d)), In view of the above discussion we can write the 
energy contribution ^ ^ due to electro-chemical changes 

as 

2 

Now Pauling's rule states that the el ect 2 x»— chemical contri- 
bution to the heat of solution is proportional to the square 
of the difference in the electro-negativities of the consti- 
tuents, On the basis of similarity of this expression with 
that obtained from density functional theory. Hodges ard 
Stott ^ have' attributed a precise meaning to the concept of 
electro-negativity by identifying it with the chonical 
potential (See Chapter l). Our attempt will be to analyse 
these data ( E ) on the basis of chemical potential. 


Ill 


The chemical potential \i can be defined for solid in 
the following way* Prom aq.n, (3.1 ) the ground state energy 
of a system B[n] can be written as a functional of its 
correct ground state density in the standard way, 

B[n] = GCn] +/ nit) dr^ + ij { 5 ) 

I r - r'l 

where G[nl = T[nJ + B„[n] 

Xw 


If we define. 


4>(r) 

^ (r) 


u M aJ 

^ 6 n- 


^ a 

^n(r) 


(4) 


(5) 


then, u = = "gr " 7 l T + 

If and the functional derivative of G are known even 
at one point inside the system then \i can be found out by 
eq.n, (5). Seitz^^ has given an expression for the uniform 
electron gas of density n placed in the field of a charge 
- £ as follows. 


= -U + - e(n) +11 ( 6 ) 

where, -U is the binding energy of the solid, is the 

chemical potential of a uniform gas, e(n) is the energy of 
a electron gas of density n and the last term gives the 
Coulomb self energy. Eqn, ( 6 ) assumes a uniform electron 
gas in a Wigner Seitz Cell. Such an assumption is 
incompatible with eqn. ( 5 ) since "if ^ (r) varies with r, then 
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2iG . 

must vaiy in a oompensating manner to give a constant 
throughout the system at theimodynamic equilibrium. It 
clearly shows that eqn* (6) is obtained by averaging over 
a ¥S cell. 

Looki n g at eqn. (5) we notice that near the edge of a 
WS cell ^(r)=- 0 and if the density here is n the chemical 
potential is given ^ 

^V-cin’nsn (7) 

This is just the chemical potential of a uniform election 
gas of density's.^. This would be the- the correcKfc p, to use 
in discussing charge transfer etc,, if the charge transfers 
were confined to a small region near the surface of the ¥S 
cell. However, our calculations show that this is not the 
case. Therefore, it makes better sense to define the electro- 
chemical potential using (5) and averaging ever the unit cell. 
For the case of an empty core pseudo potential, this gives 

p = p^(n) +1^ [ I (8) 

where H is the atomic radius, and r^ the core radius of the 
pseudo potential. The electro-chanical potential now 
depends not only on density h but also on the ionic or 
core radius r^. We plot in Fig. 4.16 p^(n), as a fiinction 
of r^, and also p calculated aoooitLing to Eq, (8) for the 
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case Z = 2, In the latter case, r^ is varied to give 

different equilibrium values of H and thus r (= 

s 

We notice that while the general trend of variation with 
Tg is the same in both cases, in the latter case the ionic 
attraction makes a large negative contribution for small r_ 
(and therefore, small E and r_). As E increases this con- 

Q 

tribution decreases (See Bq, (8)) and the chemical potential 
approaches that of the free electron gas. The large di- 
fference between the two shows that the role of the ionic 
potential is quite significant. 


In Bigs, (4.17) and(4,18), we showAQ, (p.^ - and^^NHj 
for the case of a fixed impurity (with = l) in a matrix 
(Z = 2) of varying size. In Big, 4*17> we show the charge 
transfer to the solute call (in units of e) as a function of 
The charge transfer is for the tanh case, and has 
been calculated by integrating the densiiy n(r) upto r = E^ 
(i.e, over W.S, cell where the solute is located), and 


subtracting the solute charge 
chemical potential calculated 


Z^, pg is the electro- 
according to Bq, (8)* dearly. 


the physically appropriate p^ corresponds to that confi- 
guration of p^ which has the atonic radius Eg, We notice 
that the relation AQ =^(Pg - P^) is not satisfied. How- 
ever, the charge transfer does increase with )increasing 
(Pg - p^) , and there is a large linear region. We also 
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find that rather large charge transfer (AQ 0,4) takes 
place in these model systems, Ihe reason for charge transfer 
to solute is the follo'wing. The kinetic energy term favours 
transfer of electrons from a high to low density region. 
Transfer to solute reduces the density iDhoctgnnity , and 
is favoured by the g 2 :*adient term. The electron electron 
repulsion is comparatively less if charge flows away 
from high charge density to a low charge density region. 

The ionic potential of the solvent acts in the opposite 
direction. It is not dear to us, from this model system 
study, that the charge transfer can be quantitatively 
obtained from the electro~chaiioal potential difference 
from a simple formula such as Q • 


Figure 4,18 shows log^^^ J plotted as a function 
of log^QA Q, We, remark that both of these are theoretically 
calculated quantities, and the electro— chemical potential 

B 


does not intervene. We find that a 


,4 


© # 0 < 


does go up with A Q, but not quadratically as predicted by 


simple theories and conjectures. 


In Pig. 4,19 we plot another interesting systematic 
feature, namely AQ for the case where the matrix is fixed, 
impurity or solute size is varied by varying its 
pseudo potential core radius. In this case, if we imagine 
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the solute to be oompressed or expanded to solvent atomic 
radius, the chemical potential difference would be constant 
if it were to depend only on electron density. In that 
case,AQ should be constant. Ihe fact that AQ varies widely 
indicates that the ion potential contribution to p. is 

A 

important. Since the only thing that varies is r^, and 

a (see Bq, (9)) we expect that A Q a (Pg - P^) 

a(-aA ) and would decrease as r"^ increases. This is indeed 
c ^ c 

what we find. In Pig, (4.18) we have plotted A Q r^ 

o 

for several systems, i,e, M3 - 12, M3 - H aJid M2 - H. Both 
the step model (marked 9) and continuous model (marked c) 
results are shown. The charge transfer is more in the tanh 
or continuous model because of the possibility of charge 
relaxation in attaining a continuous distribution, 

4,4 Ihe Pixiblems Encountered in Alloys with Higher Valence 
Solute i 

Pigure 4, 7(a) and (b) shows the results of fixed matrix 

and varying solute results for Ml - 12, ML ~ 13 and M2 — 13 

o 

systems where although Rg in the step model and 

r^ iC in the continuous model theAHj is negative. This 
is not acceptable because the variationally arrived at 
electron distribution shows that the matrix is intact and 
the solute cell charge density has undergone a redistribution 
there being no charge transfer, How if we 



121 


liad sliarted with, the correct ground state charge density 
of the solute then any re-distribution^eno^By of charge 
within the solute cell would lead to an increase in energy i*e* 
would be greater than and hence A Hj must be 
positive. 

To correct this situation we attenpted to start with 
a better gi^ound state charge distribution of the solute 
given by the cosine model. The alloy charge density was 
taken as 

n(?) = n° - Bj_{l + COS itr/R^) © - "? ) + 

B (l - tanh 

where is the amplitude of the cosine teim of pure solvent 

charge density in our cosine, model. This n(r) giv^ by 

eqn, (9) has three floating parameters B, r^ and ^ one of 

which can be detemined from charge neutrality and the 

other two are determined variationally , This calculation 

was conducted for Ml - 12 set and the results show that 

the continuous model Ml - 12 curve of Big. 4, 9(a) is pushed 

up (see Big. 4,20 where the two are compared) but still 

we get a negative region upto r"^ = 1,1 with the maximum 

depth at r"*^ = 0,8, This probably might mean that if the 
o 

ground state charge density is calculated in a more 
accurate model the physical inconsistency of the results 
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nay be eliiainated. 

The variational paraneti^tsiation* is perhaps quite 
unsuitable for because in these oases there is 

probably a long range charge transfer extending in the 
matrix nuch beyond the next nearest neighbour. Our 
numerical scheiaes search for a minimum in general up to 
^st ^0 ^ ^ some cases we find a minimum 

with r^^ located well inside the core region of the sol- 
vent. Unfortunately, because we are working with a pseudopoton- 
tial Such' results ar^ not reliable. 



Chapter 5 


Kohn-Sham Method in Metals and Alloys 

In this chapter we describe how the Zbhn—Shani tiethbd 
(presented in Chapter 2) can he reformulated in thb case of 
metals and alloys as a system of weakly interacting pseudo- 
atoms. In contrast to the parametric ■variational Scheme, 
here no assumption is made about the form of n(r) which is 

calculated directly by solving a Schrodinger like equation 

iq 

obtained variationally by Eohn and Sham from the density 
functional for the ground state enefgy. The nCr), thus 
determinedi defines the system uniquely and helps to find 
out the different ground state properties. 

In Section 1 we visit the concept of pseudoatom 
again as an extension to the discussion contained in 
Chapter 2, Then we present in Section 2, method for cal- 
oulating the ground state energy and hence "the cohesive 
energy of pure metals. In Section 3 we describe the 
scheme for calculating the heat of solution of alloys. The 
last' section contains the preliminaay results obtained for 
the cohesive energy of AL, Mg and the heat of solution of 
the AL-Mg alloy. 

5.1 Pseudo-atom in Metals i 

This idea as already explained in Chapter 2, is that 
a metal may be treated in some respects as if it were an 
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assembly of neutral ’atoms'. These atoms aire not the same 

as the objects we get if we take the metal to pieces 

and put the right number of valence electrons around each 

ion* In the metal, an ion is positively charged, but it 

is immetysed in the gas of conduction electrons, which tend 

to distribute themselves so as to form a screening cloud 

about each ion. The basic contention is that the object 

composed of the ion and the screening charge cgn he treated 

as the unit for building the stiucture, This^ is how 

Ziman^^ puts forward this concept in his pioneezung article 

although this idea was shared for quite some time by many 

other research workers like, Heine, Harrison, Phillips^’^, 

Prom the view-point of the density matrix, the idea has been 

44 

developed by March and co-workers , 

The above concept does not mean that the conduction 
electrons are localized* They form the usual highly degene- 
rate Permi gas. But the scattering of electrons from an 
array of pseudo-atom ordered or disordorod can be treated 
without any need to consider further shift of charge to 
maintain local neutrality. Moreover, inspite of the strong 
potential inside the core the scattering from each 'pseudo- 
atom' is quite weak and so may be treated as a pertur- 
bation due to a small 'effective potential' carried by the 
atom. If this effective potential is known one can set out 
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■fco csalculat© the band Bt27u.ctu.re, the shape of the Fe27Dii— 
surface, resistivity etc. poToperties as shown by Ziman"^^, 

Despite its simplicity this kind of model yields a 

self consistent scheme of calculation as will be shown 

in the next section. It has also been shown by the calcu— 

9 

lation of Dagens that the effective inter-atomic potential 
or pseudo-atom— pseudo— atom interaction is small as expected. 
But one of the defects of this model is that it assumes 
free electiTon states for the valence states, A nearly free 
electron fo 27 mulation of this aspect is possible but poses 
a good deal of difficulties in achieving self-consistency. 
!Ehe linearised supe 27 position that is used in const 27 Ucting 
the system has got some errors associated with it but it 
will be much smaller than that will arise when one tries to 
do screening after disturbing the lattice. 

With the above preface to the treatment of metals 
by this model we piTOceed to present the fomulation for 

cohesive energy of metals in the following sections, 

* 

5.2 Cohesive Energy of Pur© Metal t 

Here we find the cohesive energy of pure metals using 
a heutral pseudo— atom model. The basic idea is to express 
the total energy as the sum of individual pseudo— atom 
energies and their pai 27 wise inteiTactions, To make this 



I 


127 

approach practicable one has to formulate it in a precise, 
i*apidly convergent scheme. Below we present our version of 
this formulation which has been used in our study. 

Prom eq.n, (14 ) of Chapter 2 we can write the total 
energy , E[nj of a system of electron gas with density n(r) 
as 


B[n] = T[n] + +/v^^(r) n(r) dr + 


erfc 




r - r* 


ii 


( 1 ) 


where is the elect 2 ?o static interaction between the ions 

and the other terms as defined before. So, if one knows 

n(r) correctly and the functional foimis of T[n] and B [n], 

xc 

the g 2 X)und state energy B[n] is completely defined. If 
one conceives the electron density in each neutral pseudo- 
atom to consist of a uniform part n^ and a deviation fin, 
from it on account of non-uniform charge polarisation due to 
screening one can write the electron density of the system 
as 


n(?) = n^ + I fin (r - 


( 2 ) 


This is same as eqn. (28) of CJhapter 2, Here the periodic 
density in metals has been Sf)lit into a uniform part and ’an 
oscillatoiy part 6n(? - R^) centred aromad each Wigner- 
Seitz cell at lattice vector positions R^*S, Since each 
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pseudo—a'toni Is neuliral 11; follows thai; 

j an (r - R^) dr = 0 


(3) 


How the problem reduces to finding out 6n in one 
neutral pseudo— atom in the metal and thereby to find out 
the cohesive energy from eqn, (1). In both of these 
stages approximations are made. In calciaating 6n(r - 
it is assumed that only the ionic potential due to centre 
i is non— zero and the total energy is assumed to consist of 
only one centre and two centre terms. The significance of 
these assumptions will be dealt with to give a dear undeiv 
standing of this method. 


Following eqn, (24) of Chapter 2 one can write the 
effective potential as 

The approximation involved at this phase is in writing 
V[n(r)] = + I 6n(r - H^)] 

- \ + 6n(r - E^)J 

The correction term is given by 


(4) 


(5) 


Y[n(r)] - I Yin^ + 6n(r - :^)J = I , 
i " 


a ^7 




-) ( 6 ) 
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The odission of this seooncl— order terra in potential will 

affect the wave function in fourth order, i,e. 

2 A. 

0(6n^ 6nj ) or 0(6n) • This oan safely he neglected if 6n 

is snail, or if the overlap between 6n. and an. is sraall. 

^ 3 


The neglect of all other centres excepting the centre 
i can be understood easily in terns of Green's function. 

If we define 


G(r,r',w)= / < {iJJ*'(r,t) tj) (r' ,0)} > e 


-io)t 


dt 


(7) 


Then the density of electrons is given by 
^ » 

n(r) = ^ ^ G(r,r',w) dw (8) 

If we write G(r,?',w) = G°(r,r',w) when = 0, we can 


calculate G in terms G*^ as shown by the diagrams below*. 


+ ,'f\ + 




- + 


(9aj 

(9b) 


The effect of assuming all other potentials absent is 
equivalent to neglecting all the tems of class (9b) or 
multi-centre scattering tems. But the scattering fiom 
centre i is taken care of to all orders. This approximation 
will lead to an error of 0(7’^^ ) in 6n and hence an error 
of 0(7^j^) in energy. So, we oan conclude that it is 
reasonable to calculate 6n(^ - H^) in this one centre 
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appro xijnat ion . 


Following eqn^ (4) one can write 


V ^[n(r)J - I ,(r) + V [n(r)] + / 


r f'l 






r - r* 


V A. T r- ' , -*■-*■ V f ) -9- 

= { +I[7,,(a,«n(r.fip] - T^^[nJ] + H ^ ■ a^. 

1 ■ - X - - - - i |r-r' I; + ^^qL^q] - 


= I[V„,,Cr) + 6Tt„ + 1r'] + 


I'-’axi 


( 10 ) 


wliei^ n^ is the uniforta density in the i”^ Wigner-Seitz cell 


and n^(r' ) is the density of electrons in that cell, 


The tem ^ ©in. (10) can be taken care of while 


calculating B[n(r)J* It corresponds to a constant energy shift 


and so does not affect the calculation of 6n(r--R. ) which is to 


be used in solving the Schro dinger-like equation expressed 


by eqn, (25 ) of Chapter 2, 


h ^(r) = e ^^(^ 


( 11 ) 


where 




[nj + 6n(r-fij^)] 


(12) 


r - r* 
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Now 6ii(r - is given "by 


6n(r - fi^) = ^ I ^\r) | - n 

Xocc 


(13) 


This three equations iron (ll) to (13) can he solved self- 

■> 

consistently to get 6n(r - E^). 

Por should have the ionic potential, i,e,, 

suppose for the appropriate ionic configuration, the wave- 
function hence the charge density known. Then, 

the coulonh potential is given by 


=f 

* coul ■' 


(r* ) 


ion 


■*‘1 1 

r - r' 


dr' 


(14) 


In addition, the conduction electivons and the ion-core ele- 
ctrons being both identical femions will exchange and oorrs' 
late. This contribution is given by 


" 6n(r-^)] (15) 

Now the ^Q 2 t: gi'ven by 


7 ^ ,= 

ext coul xc 


(16) 


- But in our study we do not calculate the 
procedure. Instead we use a pseudo potential whose choice 
will be described in the next section. This pseudo -potential 
simulates the effect of eqn, (l6) in a given 
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ana E[nJ = 

nn dr dr» 


- i/ 


o 0 , 


T- + 7^^ 


Cl9) 


f r - r^ 1 

Subtrac1.ijag: eqrii (19) fron (18) we will be left with the 
AUj'b and ATJjj's. First let us evaluate the change in one 
centre teims ^ phase shift analysis we can write 

T (e^ - e°^^- i J f ^ j 


-ft ^ i I('‘^^+l)C£t(ep)-6^(0)]ej 


0 Ep 

” it I ^ ^ +l)5^(e)de - ^ 


( 20 ) 


In the case of pure netals A Z = 0, so only the first tem in 
eqn, (20) survives* 

ITow we Set out to find the change in the exchange and 
correlation terms given by 


I(^^Q[n(r)] - •M- 2 .QCn(r)J} n(r)dr -/{e^Jn^j - ] } n^dr 

Sorting out the one centre terms from this expression we get 
a contribution of 

^ / (6ei^ - 6^^) [n„+en(?-^)]d? 

where = e^Jn^.h5n(^-^)] - e^^[n^J 

and 6^1^^ = n^^[n^.rtn(;-^)] - 


(21) 
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Tie change in coulomb tem can be written a^ 

- i ! 4^ 4?' + i if aJ. 

I !■ - r'l . I? - ?'| 

and after sorting out the one centre terms we get 

6n^(r) 6n^Cr') 
Ir - r’l ' 


.In/ 

n* O I 

1 r - r’ [ 


dr dr’ - I dr d?» 


We write these terms in Fourier space as follows i 

6n(r - iL ) = ^ 6n§ and = I -iS- |r-r' | 

1 ^ \x^P\ \ 

where is the volume of the space containing r and r ’ f 
the first term in the expression (22) becomes 

-% L H d?d^' 

<1 5’ <1^ ^ 


-n_ I 


47c 6n|^ 


/eil»r 


° 4 4" 

6n^ 

= -4u n ( — a) 

q, 0 q.^ 


dr 


27^ ^ 

= + — / r^6a(r) dr 

The Second term of the same expression is 

n I 


( 23 ) 


( 23 ) 


( 24 ) 
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Now we can write the expression for the one centre teim as 


‘% \ (2^+1)/ ^6| (e)de + /(6e^^-6ti^^)[nQ+6n(r)]dr 


2 %n 

■^3 


o f 2c V (6^?r6n-q) 

— /r^6n(r)dr - fi I — 1- 


(25) 


1 q.* 

The two centre terns A 15:1 can be estimated from the fun- 
ctional relation of kinetic# exchange and correlation 
energies 

T[u^+4n(?)] + E^Jn^ + 4n(?)] = + 


■ an 


/ A n(r)dr 


n=n. 


+ i / An(r) A n(r ' ) 






9n(r)3n(r') n=n 


dr dr' 


(26) 


As shown by Eohn and Sham^ one can write 

3 (l+\c^ /-► -► N 1 T -i "?•(?-?') 

— ■ I = E(r-r') = ^ e 

9n(r)9 n(r') n(r)=n^ , 1 1 

where K is the response function and ~ is the 

dielectric function and An(r) = I is a peiiodic 

i 

function having the lattice periodicity* 

Pirom eqns* (18) and (26) we can collect the two centre 
terns to find outAIJ^j as 

^ Ujj = J 6n(r-fi^)6n(r-B.^ )Z(r-r’> dr dr' 

6n(2>-Bj^) 6n('r'-$^) ' 


+ i / 


r - r' 


dr dr' 


(27) 
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.f 


tfri'bxng 6ii and K in terms of Fourier series expansion 


6n(r-:^) = I 6nJ 


and 


X(r-r') = Y E e 

i 1 


-iq, (r-r ' ) 


¥e get. 


AH = ilKJ 

q 2 ;^2 et-1^ ® ^ ^ 


q‘ "q^ 


( 28 ) 


Tile oont ribution of tiiis term to tlie ground state energy 
can be written following eqn. (17), as 

^ ^ ^II ~ ^ 

z 


where = (^ (| \-ilp = I 

91 9 . 


(-iK )^iq.|R Ejl 


This can be looked upon as the interatomic potential between 
the ps eudo~*atome at lattipe ^sites i and j,' This expression 
unlike its counterpart in the pseudo -potential theory, does 
not use the linear response approximation, i,e. the assump- 
tion that the electron gas respoitds linearly to the potential 
The treatment is exact apart from the local density 
approximation for the exchange and Correlation contribution 
to the potential, assumption of a spherically symmetric 
potential, and the neglect of terms in density arising from 
two centre interference, and neglect of three and higher 

centre contributions to E , 

S 
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In the treatment of Dagens^, a neutralizing auziliazy 
charge distribution n^(?) at each lattice site is postu- 
lated, (This has a Q-aussian or similar form, anri = 2 • 

The effect of the difference between and the unifom 

distribution n^ spread over the Wigner-Seitz cell is treated 
in perturbation. Assumption of a neutral auxiliaiy pseudo- 
atom is unnecessaiy ■ as shown above; it is a source of 
arbitrariness since the distribution n.(r) is not unique. 

Ball has used n^ over the ¥S call directly, as we have, but 
costs his theory in the form of pseudo -potential perturba- 
tion theory. We have shown above how the density functional 
formalism can be directly used, once the pseudo— atom charge 
density is known, to calculate the energy change due to the 
^ coupling of pseudo-atoms via the electron gas. Since the 
6n(r-R )'s are localized around their respective cells and 

" A ^ 

Since the coupling kernel E(r-r') is short mnged 

(1 r-r* ), the two centre term can be expected to be 

' * 

Small, 

The advantage of this method is that it is not necessary 
to use a pseudo-potential-. . The two centre or inter- 
atomic term being small i,e, of the order of 4 percent of 
one centre term the method is rapidly convergent. On the 
question of whether- it can be app3J.ed to transition metals 
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or not good prospects appear to "be there for the following 
reasons: (a) The transition metals also form close packed, 
non—directionally bonded structures, (b) The eneigy differen- 
ce between bcc, fee and hep structures are small confirming 
non-directional bonding (c) The inter-atomic potentials 
derived from the structure factor data appear to be q.uite 
strong in transition metals indicating stronger intercell 
coupling but again it might still be small as the ratio 
of this to the cohesive energy may be as small as in simple 
metals(transition metal cohesive energies are much larger 
than in simple metals the ratio of melting temperature to 
cohesive energy ) seems to be the same for say 

the alkali metals, and noble and transition metals). 

The identification of the two centre term with the 
interatomic potential implies physically that each ion 
caxries its neutralising Wigner-Seitz cell of electrons 
around it. While this is acceptable when the ions are in 
lattice positions, in a liquid or vibrating solid this premi- 
se has to be examined. If this is correct there also then 
we have a simple method for calculating the interatomic 

potentials which appears to be valid for transition metals 
as well . 

Using a spin density functional one can try to find the 
condition for the pseudo-atom' to be magnetic. This also will 
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give the spin density oscillations, hyperfine fields 
etc. Ihis area is still totally unexplored through this 
method. 

5.3 Ihe Heat of Solution of Alloys: Simple Alloys : 

Suppose in a system one of the lattice, site i has ion 
A with charge whereas all the other sites are occupied by 
ions B with charge Z^. Ihe energy of the system can be 
written as 

B = + (H-l) Eg + (30) 

( B ) 

where E^ ^ is the energy of the pseudoatom A as it exists 
in the alloy. Eg is the energy of a single pseudoatom B in 
pure B and B^ is the interaction energy between pseudo atom 
A and the surrounding B ps eudo atoms if similaiy 
define B^ aS the energy of a single pseudo atom in pure A, 
the heat of .solution per atom of A can be written as 

A H = B - E^ - (H-1) Eg = b|®^ - ^A ^^AB 

How let us tiy to write the total energy of the system in 
the form of eqn. (30). Following eqn. (17) if we wiite the 
total energy as 

B = ECn^] + A Ug + AUgg (32) 

with n^ as the average density of electrons in pure B and 
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A TJj and A Ujj as one centre and two centre terms respectively^ 


we can express E[n^] as 


E[n^J ={l[nj + + 

(H-l) i I 


l5^ 0 I % I 


H-l 1 ^ 


Zb A a 


B l3^ 0 I Rj_ I 

where the superscripts like N-l, 1 indicates the number of 
pseudo atoms, the subscripts identify the cells, is the 
electrostatic enei^ of pseudoatom B as in the pure state* 
Similarly we can writeAUj as 




JSB 1 




O'eB 


■J . 


+ I / 


jsB 

-il II 

jeB 


f r - r» I 

6n®(J-R.)6 ^ ^ 

^ dr dr' 

1 r - r' I 


dr dr' 




t| I (21+1) 


de 


de 






= (U-l) AU® 4- A 


(33) 


where fin is the deviation of charge density from the uni- 
foim for pure B and fin^ is the deviation from n^ In the 
Wigner-Seitz cell of solute in the alloy. Ihe two centre 
tern also can be evaluated in the same way to give 




j 


B V . 

yi 6n^;)finB(?.-5.)t^^ ^ 2rr-?.) Jd? d?- t ^ (34) 


The first Wo teims of eqn. (34) can be writt 


en as 


)dr - 2 ^ '7®(r)6n®(r-R. 

+ I /an®(r-^)6nBCJ^H )[ — ^ + K(?-J')] d5 dj' 

J 0^ + E(?-?')]aJ d?' 


(34a) 
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This can be better understood in teims of bonds because the 
difference in the number of B-B bonds if one site is 
occupied by A with that if it had been occupied by B is 
m = I Bo, if IS ths bond 

energy per atom in pure B the expression (34a) can be 
written as 

H bBB - 2 B®® = (H-Dllgj - Bj3 

The rest of the teims in eqn, (34) can be written as 

A 

T? - y f d?»3d? + 

|r - r'l 

T r r 611 (r) ^ 

+ I C / - — # dr J 

- jj^.O i ^ 


A Z Z 

+ 2 I rdn^(r) E(r,-r') dr dr' + " 3 (34b) 

- i-^9 ^ ■ ■ I I 

The evemation of this type of terms have been discussed in 
the pure case and the only difference here is J6n (r)dr = 
Z^ - Zg = A 2* The contribution from the two centre terms 
can be written in terms of eqn. 34(a) and (b) as 


^ ^AB ” ®AB “ ®BB 


(35) 


So collecting all the one centre, two centre terms together 
we can write the final expression for the heat of solution 
per atom of solute as 
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AH = § I (21+1)/ ee ( 


P d6r(e) 

f ^ \ 




It 


0 


de 


XC ^XC' ^ 0 


6n^(r) n©(r'iig.) ^ r .6n^(?l 6n^(r' ) .p 

J dr dr» - -i- J . -^ -►- • 


r - r' 


r - r' 


+ ^ab + + \oC“o3 +K\- \ 


(36) 


The calculation of one centre tenas involved in the eqn. (36) 
can le done in the Fourier space as detailed in the disccu- 
ssion. for the pure case. 


5.4 Oalculation and Results 

Tn the present section we will descrihe the niuaerical 
methods • followed to solve the Schrodinger equation like 
Kohn-Sham eqn. (11) and present the result for the phase 
shift and 6n(r) in 41, Mg and Al-Mg alloy. These results 
are non-self consistent and the prohlenf^ of self-consistency 
will be discussed in the context of these reSTilts, From the 
6 *s and the 6n(r)’s we evaluate the cohesive energies for 
Al^Mg and the heat of solution of the Al-Mg alloy by eqns, 

(17) j (33) respectively. These results will also be presented 
here. 

(i) The Solution of the ES-equation : 

The ES equation poses an initial value problem for a 
linear second order differential equation and there are 
several methods to solve it. But to follow any numerical 
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method for its solution the first task is to divide real 
space into a mesh. Then the eq.uation is solved at each 
point of thie mesh, 

For the case of atoms, the radial wave fn is oscilla- 
tory at small and medium values of r, and exponent t a lly 

A 

damped at large values of r. Our case is quite different* 

In the begi nni ng (near origin) we have an exponentially- 
damped radial wave function and at large r it is free ele- 
ctron like. The charge density (rather the difference in 
charge density from the initial unifonn distribution i,e, 

6n(r) has oscillations through out. The determination of 
phase shifts requires accurate values for the radial -wave 
f-unotion at large r where only can it be well approximated 
by its asymptotic fom, 

From these considerations we can see that our choice 
need be different from that of Herman and Skillman^ ' 

for the atomic case. Our mesh is uniformly spaced at distances 
of A r and the KS equation is solved upto a distance of 
J_ ■, of where the phase shift is determined by matching 

it with the asymptotic form. In general, smaller r and 

larger give better results. But the limitation of the 

memory of the computer demands an optimum choice. 
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The method followed to solve the ES equation numerica- 
lly is due to UumeiHDV which is descrihed in det^l in the 

49 ■ 

book by Hartree -r, /• 

•It runs as follows: for small r(i,e. 

the first two points on the mesh) the KS equation is solved 
analytically for a given potential and a given energy. 

Since the potential in the cases considered is highly posi- 
tive near origin the solution will be in terms of modified 
spherical Bassel function (we use the condition that the 
radial wave function is finite at origin), After knowing 
the solution at the first two points of the mesh we start 
the outward integration procedure at each outwar^i point 
involving the adjacent iwo inward points where the solution 
is already known (see the book by Hartree ), This was 
done by a modified version of the computer programme written 
by Hermann and Skillman for the case of atoms. 


Since we need the phase shifts for our purqjose we had 
to incorporate this calculation in the above computer 
programme. The asymptotic form of a phase shifted radial 
wave function for the angular quantum no.l is given by 


Lin Ej_jj.(r)5=cos 6^ j^Ckr) - sin 6^n^(kr) (37) 

r large 

where and n^^ are the spherical Bessel and Heiamann fun- 
ctions for a given 1 and the energy B = k^ in the 
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Hydlergs. If we choose test radius at sufficiently 
large r where the radial wave function has taken the foim 
given by eq.n. (37). Ihe phase shift can be expressed as 


tan 


kn[(kB^) - 


(38) 


where 

T 

Y = ~ ^ can be found out firom the numerical solution des- 

'1 H]_ dr 

cribed in the last paragraph, and n^ are derivatives of 

3^ and with respect to r. 

To choose an optimum mesh size r and test radius 
for a given potential we vary both these quantities ^ r 
and and study their effect on 6 *s, it r was varied from 
.01 to .08 a.u. and B^ was varied from B^ = 4.0 a.u. to 
14,0 a.u. The 6£Swere found to be stable for^r = .05 a.u. 
and B^ =6.0 a.u. So, we choose for our calculations a mesh 
size A I* = .04 and deteiminie the phase shifts at B^, = 8,0 a.u. 

After calculating the phase shifts we need to calculate 
the deviation 6n(r) from the starting distribution of the 
electrons (which we assume to be uniformly n), 6n(r) is 

given by 
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where the energy E = and ICp is the Perrai nomentum. Ey 

e^ps^hms n in tenas of spherical waves and writing 

®^licitly in tenas of the radial wave functions 
one can rewrite the eqa, (39) as 

Tf CO 

6n(?) = ^ p k“ dk I (21+1) { [E,^(r)]= - jMkr) > («) 

0 1=0 . . - 1 

Once values are calculated over the mesh, at each point 

the 6n(r) can he determined hy eq.n« (4. The normalisation 
of values obtained numerically can be done by comparing 
it with the etin* (37) which is given in the normalised foim* 

With this 6n(r) given by eqn. (40) one can find the 
modified n(r) (let us denote it by n 2 ^(r) of the system at 
the end of first cycle as given below 

n^(r) = n + 6n(r) 

This n^(r) can how be put in eqn. (16) to give the effe- 
ctive potential to be used in the next cycle to solve the 
S3 equation. This way one can go on until the self-con- 
sistency is achieved. The final values of the 6j^(5)'s for 

->■ ' 

the given energy s and 6n(r) obtained self— consistently 
can be used to find out the cohesive energy of the pure 
metals and the heat of solution of alloys from eqns, (17) 
and (33) respectively. 


The above calculation should give self-consistent 
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5 2^3 whioh. satisfy tlie Priedel sum rule. 

^Z. = I i (21+1) 6-,(e^) (41) 

^ 1=0 ^ ^ 

where e^, is the Eermi energy of the system, AZ = 0 for 
pure metals and equal to Z^ Z^ for the alloys, we can 
start from any form of initial n(r), the self consistent 
solution will give the same final n(r), and satisfy eqn. 
(4l)» But a better choice of initial n(r) used in calcu- 
lating will converge to the self consistent solution 

faster. Since we find from our experience with the parame- 
tric variational scheme that the uniform dhaiige distribution 
gives a good approximation for the ground state charge 
density. Our initial charge density is taken to be uniform 
(H) and to have the equilibrium value, i.e. to contain 
charge Z^ in a volume of rarius But with the Ashcroft 

empty core pseudopotential and this choice of initial n(r) 
the Briedel sum rule is not satisfied. So, we vary the 
atomic radius (which is equivalent to changing n) until 
we can satisfy the Briedel sum rule given by eqn. (4l)* 

It is quite a standard procedure to start with an initial 
distribution which satisfy the Priedel Sum rule (Lang and 
Kohn^, Popovic and Stott^^). 

To check the accuracy of our numerical scheme we find 
out the phase shifts for different 1 in absence of any 
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potential (l.e. = 0, ). 

sixth, deoimal place for = 4,0 

mal place for = 12,0 with Ar 
cases. 


This yields zero 6 ^ upto 
a.u, and upto eighth deoi- 
= 0.04 a,u, in both the 


The 1 -values are generally taken upto 1 = 2 which we 
also follow in our calculation. But to check the a^nount of 
error involved we have calculated the 1 = 3 contribution to 
the Priedel Sum rule given by eqn.. (6) as (aZ)^= 14 6 (ej,)A. 
It is found to be quite small ( 0 . 0046 ) for a calculation' 
with the empty core pseudo potential for magnesium. (A 
is also found to decrease with decreasing uniform density n,. 
The empty core model potential tends to pile up charge in 
the region r^< r<R^ and this tendency becomes more promi- 
nent with iteration. This defect of the model potential can 
be overcome by choosing a flatter, smoother potential such 
as a local version of Animalu-Abarenkov Heine potential, 

With all the precautions mentioned above when the 
calculations were carried out taking an empty core model 
pseudopotential for the ion it was found that after one or 
two cycles of relatively moderate deviation the 6n(r) starts 
changing violently. Animalu-Abarenkov-Heine (AAH) model 
potential was used (i,e, 1=0 component is only taken) 
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in its place. The local AAH potential is (as described 
in Chapter 1) given by 

^ion^^^ = “ ^ ^ 

7 / . 

- ixjx for r > r 

o 

The values of r^ taken from Shaw*s calculation,. For 
Mg 1*0 = 2,58 and for A1 r. = 2.15. 

The results Porr^6n(r) at the end of first cycle for 

Al, Mg and Al-Mg alloy are plotted in Pig. 5.1(a), 5.2(a) 

and 5.3(a) respectively. These also shows V(r) in Pydberes 

in the accompanying figures marked (b) for each of these 

r^6n(r) plots. The r^6n(r) curves for the pure elements show 

.23 

that Some charges have either flown out /for Al) or 
moved, inside (.016 for Mg) the sphere of radius This 

may be due to the fact that the phase shifts are inaccu 2 ?ately 
calculated while satisfying the Priedel Sum mle. As a 
result we might expect that the estimate of the cohesive 
energy or the heat of solution fjcom these non-self- 
consistent calculations will be in error. But still we go 
ahead with the calculations so as to convince ourselves about 
the practicability of the schemes proposed. 

The complete procedure involves finding Eg as a fun- 
ction of n^ and locating the minimum. That is, for a 





Vcff(r) in Rydberg 




eft (f ) Rydberg 
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given (and hence lattice constant) we obtain the self 

consistent 6n(r) and the corresponding Eg. This is done 

for several values of n i the lowest value of B is 

o g 

cohesive energy and the corresponding n^ is to be compared 
with the observed density. Further, in principle, the ion 
potential may be used directly (as discussed in, Section 5.2) | 
it need not be modelled by a pseudopotential. Unfortunately, 
since we have not solved the self- consistency? problem, we 
have not done the complete cohesive energy calculation as 
above* For an value given by the experiment, we vary 
the initial potential as described above till the Friedel 
Sum Rule is satisfied. For this neutral configuration (not 
self (Consistent yet) we calculate E . The calculated 

s 

value is expected to be close to the observed since 

most of the contribution is from n . 

o 

Table 5.1 shows the phase shifts 6^ 's for 1=0, 1 and 

jr 

2 for Mg, A1 and Al-Mg alloy . In the pure A1 or Mg metal we 
find that although the phase shifts at satisfy the 
Friedel Sum Rule approximately still the r^6n(r) curves 
plotted in Fig. 5.1(a) and 5.1(b) show an amount of charge 
flow ( to or away from the origin. This can only result 
from an inaccurate calculation of phase-shift which involves 
derivative of R- and hence is quite susceptible to error. 



TABLE 5*1 

Variation of Phase Shifts with Ehergy for Al.M^ 

and Al-Mg 
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Table 5.2 shows the eatimated. value of the different 
terms involved in the calculation of cohesive energy of 
pure metals by eq.ns, (24) and (25). Some of the terms which 
give positive values for are found negative in the case 
of Al, This may be due to the fact that r^ 6 n(r) has positive 
deviation from zero in case of Mg while it is n^ative for 
Al, The comparison of the observed cohesive energy with 
the calculated one is good in view of the fact that Priedel 
Sum fttle has net been satisfied. But the difference in the 
calculated and the observed cohesive energies of Aluminium 
is quite large. This may be due to again the flow of 
electrons away from the region of r < We see also that 

the charge oscillations are fairly prominent. We note that 
the contribution of n^ or the uniform part to the ground 
state energy is much larger than the residual one centre, 
and the two centre terms. The ratio is 170:3^0.5 for 
M and 38l!lJ0,2, Thus, as seen earlier, most of the co- 

s 

hesive energy is due to the uniform part of the electron gas 
density. The smallness of U f justifies our neglect of 
three and higher centre terms. The cohesive energies are 
close to the observed values. 

Table 5.3 shows the calculation of the heat of solution 
for the Al-Mg alloy by eqn. ( 35 ). 
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TABLE 5.2 

Values of Different Tezma in Equation (25) and '(26) for the 

Cohesive Energy Calculation 


Bqn. (25) J 




A1 

1st term : 

-.008244 

.093452 

2nd tena ; 

+.018768 

+.02&366 

5rd teim : 

.017928 

-.019492 

4th teua s 

-.000482 

-.001156 


.028170 

.101171 

Bqn. (26) : 

.005414 

.021417 


.053584 

.122588 


-1.704676 

-3.960880 

Calculated cohesive energy 

-1.671092 

-3.838292 

Observed cohesive energy 

-1.784 

-4.14 


All the energy teims are expressed in zydbergs 
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Th .0 first term is tiie cjiaxige in tlie one electron energy, 
i® positive^ is tiie ezciiange and correlation temi* 

The third term, representing the electrostatic interaction 
betweeA the uniform electron distribution in the matrix 
with the charge deviation dn(r) centred at the impurity 
is negative. This is because the 6n(r) is negative (there is 
one less electrons in the impurity cell centred pseudoatom, 
than in the host pseudoatom). The next tern, being the 
negative of half the e.s, self energy of the distribution 
dn(r}, is also • negative. The fifth term represents the 
electron gas mediated interaction between Mg and A1 pseudo- 
atoms. This is, interestingly, attractive. To these, we 
add the uniform -A1 contribution and subtraict the g.s.e. 
of Mg, The value we obtain is negative, but about three 
times larger than the observed value. 

We cannot draw too many conclusions from this incom- 
plete and isolated result. There is every reason to 
believe, however, that if the self consistency problem is 
Solved, the above Kbhn Sham p 2 X)cedure would lead to a veiy 
good quantitative and qualitative understanding of cohesive 
energy of metals, and heat of solution of alloys. The 
method is free from the inadequacies of the variational 
procedure of Chapters 3 and 4. It req\h.res no pseudopoten- 
tial, no assumption of form for n(f). The approximations 



TABLE 5.5 


Oalculation of, tfce Heat of Solution of Al-M Alloy 

S 

by eqn. (36) 


Bq, 33 ., ( 56 ) 1 1 st term : .8975 

2 nd tem J .7372 

3 rd term • i - 1,2140 

4tii tem : -,2321 

5tli teim ; -.0783 

6 th. to 8 th. teim i — 1,9567 

9th term -I. 67 II 


- .1753 

54.8 K cal/gm, atom 


illl the energy terms are in ^dberys 
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made seem (both, on the basis of our oal collations and self 

9 

oonsistent oaloulations of Catena for alkali metals) 
to introduce negligible errors, and to be -well justified. 

The various terms oontributi 23 g to can be physically inters 

s 

preted, and the charge distribution is also obtained. We 
expect to puriitie this approach to a conclusion. 



Oiiapter 6 


Oonoiusion 

Ve have described in this thesis an attempt to study 

properties of dilute simple metallic alloys using the density 
* 

functional method* Xhe properties of interest are heat of 
solution, cd&ajRge transfer etc. Ve have shown how the density 
functional formalism can be adapted to the case of dilute 
alloy, i,e, a metal atom A embedded in a matiaz B. The two 
Qommonly used versions of the densitT" functional formalism 
ai*e variational parametrization of the density, and use of 
the .variational principle to obtain a Schrodinger like 
equation from this the wave functions ^ , and the density 
n(,r} are calculated. The latter determines the potential, 
and there is a self-consistency problem to be solved. Once 
the density is known, the ground state energy B and iiience 

s 

the heat of solution can be obtained by using a well known 
functional form for B as a functional of nC^). 

s 

In recent years, both methods have been applied with 
great success to the pDroblem of meta3J,ic surfaces. Clearly, 
the variational parametrization method, though computationa- 
lly simpler, is limited by the choise of variational fun- 
ctions possible with a few parameters. The latter or Bohn 
Sham method gives more realistic n(r), and is much more 
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powerful. It has recently been applied to the ohasisorption 
problem (oalculation of chemisorption energy of a chaise 
Z as a function of distance d from the surface) by lang 
and co-workerS^® again with great success. 

Initially* the parametric variational method was used 
by us oh aeoount of its simplicity. The results obtained 
oan be viewed in two different ways; either in terms of 
how well they relate to real metals, or as studies of 
j^'^eiabtic trends in stable model metallic systems similar 
lb real metals. The results obtained by us do not always 
describe with reliability and quantitative accuracy* the 
observed heats of solution of real metals. There are many 
good reasons for this. The most important is the fact that 
the heat of solution is a very small quantity, typically 
^1 percent of the cohesive energy. Thus err*oirs which do not 
matter in cohesive energy calculations can be important for 
heat of Solution, For example, the periodic component of 
the density, though large, contributes very little to the 
pure metal cohesive energy. However, its presence may aid 
or hinder the process of charge adjustment in the alloying 
region, and may thus influence the heat of solution. It 
is difficult to include this effect in the parametric 
variational method without complicating it enormously. In 
the one place we have tried to include it (Ml - 12 system 
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See Chapter 4 Seotior 4), the effect is sizeable, and in the 
right direction* We have also had to use a pseudopotential 
to represent the ion {this limits the applicability of the 
method to sickle noatrjmsition metals), and have used a 
local density fuactional formalism with a gradient term* 

Por systems with 2^^ < 2-g, we find alloying with 

charg# transfer* varying the core radius of the pseudo- 

potential, host or impurity equilibrium size can be 

as a function of the host 

^ fOie heat of solution has been studiedZor impurity 
r^fcdiae, and systematic trends have been found (these are 
discussed in Chapter 4)* We have also attempted to analyze 
these trends in terms of concepts commonly used in alloy 
theory, namely size and electrochemical factors* lEie first 
problem is of quantitative definition of these. We have 
attempted some definition* However, though the latter can be 
identified with the electrochemical potential p which has 
a precise definition, its accurate calculation has not been 
possible. Two estimates are presented. Using these, we 
have tried to corrdLat e charge transfer AQ and heat of solu- 
tion (except for the size factor term) with - p^. 

Though there is a definite norrodatron, it is not as simple 
as suggested by * small deviation theories. In addition, 

* Jt 

to the uncertainty in p , a contributory cause could be 
the fact that charge transfer takes place well inside the 
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impurity VS oell, and not on its surface or over a large 
volume, IHius, local effe^sfes may be significant and over- 
all (Quantities like \i may not be tbe only important factors, 

Ve b-ave al^o presented, inter alia, calculations of 
for tb® Si pure simple metal in the local density 

approximation in wbicb tbe periodic component is not calcu- 
lated in perturbation tlieoiy* These oal cxilations show the 
periodic density component to be substantial but its contri- 
tnition to g.s.e, to be small. 

In Chapter 5» we describe how the Kbhn Sham method can 
be reformulated as a weakly interaoting pseudoatom model for 
pure metals and dilute alloys. As mentioned earlier, in 
this method, no assumption is made about the fo23]i of n(r) 
which is obtained from solutions of a Schro dinger like equa- 
tion, This method, as formulated here, (can be expected to 
lead to reliable values for alloy parameters, as well as 
to qualitative understanding of the values. We have not been 
able to solve the self-consistency problem here, and so 
present only preliminary results of a single calculation, 
for an Al-Mg alloy. 
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The linearized integral equation obeyed by the irredudble vertex function associated with the density 
fluctuations is solved exactly up to second power in the wave vector. This is used to compute the static 
polarizability of the h(»nogeneous electron system up to this order. It determines the gradient expansion 
coeffiemt in the density-functional formalism for the inhomogeneous electron system. The result is compared 
with existing approximate caclulations. The method of solution is applicable to irreducible vertex ftmetions 
which appear in the determination of other oorrclation functions of Ihe homogeneous systems. Our gradient 
term vanishes for both extreme high- and low-density regions, unlflce the results of Kleinman and Sham. 


I. INTRODUCTION 

Almost a decade ago, one of us^ developed meth- 
ods of investigating the integral equation obeyed by 
an irreducible vertex function associated 

with the transverse spin susceptibility of an itiner- 
ant magnetic electron system. In the long- wave- 
length limit q-*0, it was found that one may expand 
this vertex function in a power series in the wave 
vector, the coefficients of each power of which is 
a spherical harmonic expansion. In this way, the 

coefficient of the spin- wave dispersion is found 
to be an integral over where the super- 

script Indicates the coefficients of ^ in the 
expansion and the subscripts 10 mean the ^-wave 
part with zero projection in the spherical harmonic 
expansion. It was then found that obeys* a 

linear integral equation which in turn was solved 
by a variational method, since a series solution in 
k of this equation was not found satisfactory. In 
this manner, the effects of interactions (including 
all renormalizations) among the electrons are 
taken into account almost completely within the 
scheme of the linearized-vertex-function formal- 
ism. To investigate the spin-density-wave (SDW) 
instability of this system, the solution of the static 
part of this same vertex equation for finite wave 
vectors is required. To examine this, we developed 
a variational method for T^,[k;q) directly and de- 
duced from this a simple criterion for the SDW in- 
stability. The latter method was found very power- 
ful in dealing with other vertex functions, and the 
author with his collaborators has examined other 
properties of the interacting electron system, 
as well as its relationship with other methods such 
as moment- conserving schemes.*'® 

In the density- functional approach to the inhomo- 
geneous electron systems, one finds that the bet- 
ter the knowledge of various correlation functions,® 
the better the description of the inhomogeneous sys- 


tem. For instance, in Ref. 9, we used the Yukawa 
interaction model in the limit of large screening to 
deduce the spin splitting of energy bands in ferro- 
magnets. In the paramagnetic case, in particular, 
Sham^° has given a calculation of the static density- 
correlation function up to terms by solving the 
integral equation obeyed by the associated irreduc- 
ible vertex function by one iteration. The purpose 
of this paper is to use the methods of Ref. 1 to 
determine the same quantity that Sham computed 
both by an exact solution of the vertex equation and 
by the variational method for which the general 
expression for all q for the quantity of interest al- 
ready exists.® In contrast to the spin-wave calcu- 
lation,. it is found here that the vertex equation can 
be solved exactly in this limit. All the integrals 
appearing in the theory can be done generally for 
a wide class of static interaction potentials. The 
detailed calculation for the Yukawa potential is de- 
rived from the general expressions. This there- 
fore serves to determine how good the simple vari- 
ational answer is to this order, as well as derive 
an exact expression for the coefficient of the square 
of the gradient of the density in the theory of the 
inhomogeneous electron system. KLeinman^^ has 
very recently developed a solution to the vertex 
equation to this order. He has criticized Sham’s^ 
work in the light of the results he has derived. We 
will here make a critical assessment of the situa- 
tion. 

In Sec. n, we give only a brief description of 
the solution to the vertex equation appropriate to 
the problem at hand. Elsewhere^ we will give a 
complete account of this method and its application 
to other properties of the system. In Sec. HI, the 
variational solution® is given to this order. The 
results are compared with each other in Sec. IV. 

In Sec. V, we specialize the results for the Yukawa 
interaction model, and we give a critical assess- 
ment of these expressions in relation to the works 
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FIG. 1. Graph of fi as a functioD of irkp in Hartree 
units; Sham's: dashed curve; Klelnman’s: dot-dashed 
curve; exact, Thomas-Fermi: solid -triangle curve; 
exact, self-consistent: solid-line curve. 


of Sham^° and Kleinman. In this section, we also 
point out the adequacy of the schemes for calculat- 
ing the screening parameter in the Yukawa model. 
This is because with a pure Coulomb potential for 
the electron interactions, there is a divergence 
characteristic of the long-range nature of this po- 
tential. A standard way to handle this is to 
introduce a screened Coulomb potential instead; if 
the screening is treated as a density -independent 
parameter, in the limit of it becoming zero, this 
divergence persists as Kleinman has shown, which 
is also borne out by our analysis. The forms of 
the screening as a function of density, however, can 
be treated in some schemes which are discussed 
in Sec. V. The quantity of interest in the gradient 
expansion is Q (in the notation of Sham^°) as a func- 
tion of density, and it is given in Fig. 1. For com- 
parison, we have also plotted on the same graph 
those given by Sham and Kleinman. As shown in 
this paper. Sham's analysis is exact for high densi- 
ties and is zero in this limit; we are in agreement 
with this result in contrast to Kleinman, For low 
densities, Kleinman's analysis is correct within 
the density -dependent screening model where Q, 
again approaches zero. Sham's procedure is in- 
applicable here, as his analysis was specially de- 
signed for high densities, as shown in Sec. V. 

These features are brought out explicitly in our 
analysis. Section V also contains details of these 
calculations. These various features are due to 
a proper combination of the renormalizations and 


to the screening being a function of density in the 
form specified in Sec. V. Even more important, 
the antiparallel spin correlations are not taken into 
account in a linearized vertex equation, and only 
recently have we made progress in a variational 
solution to the nonlinear vertex equation. We 
also present here a short summary of the results 
obtained. 

n. SOLUTION OF THE IRREDUCIBLE VERTEX 
EQUA.T10N 


The longitudixxal static dielectric function may 
be expressed in terms of the irrcxiuciblo density- 
correlation function and the associated irreducible 
vertex function as 

€x(f7) = l+ (4n/^^)x(^) » (1) 

where 

x(q) = 2e^ f JF(fe, <?)r(fe;^)[rf»fe/(27r)*l • ( 2 ) 


Here £F(fe, ?) is given by 


£F(fe, ?) = 


/o(^+ ^q) aq) 


(S) 


with the noninteracting one- electron energy, attd 
/o(*) the Fermi function associated with the inter- 
acting electron of state It. r(A‘;(7) obeys the linear 
integral equation in the random-phase approxima- 
tion when the antiparallel spin correlations are 
neglected, '* 




This equation is the same one that Sham*® and 
Kleinman** considered except that we have written 
it in such a way as to display in an explicit way 
the correlation contribution to the one-electron 
energy. It is for this reason we use the notation 
r instead of their A. The advantages of doing this 
willbecome evident shortly. The interaction poten- 
tial "U ,( I E ~ Ej I ) represents the statically screened 
electron-electron Interaction. It Is for this reason 
r(fe;< 7 ) does not depend on the frequency part asso- 
ciated with k. We shall not specify its form here. 

We will now express these equations in dimen- 
sionless form, as is done in Ref. 1. Thus, in the 
notation of Ref, 1, we have 


x(<r) = I 

and 

r(x;.v) = l + ^^^^| Y,(|x-Xi|) 
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Once again, provided that Eq. (14) holds, we obtain an explicit solution for Vf^hc) which is nonzero only 
for Z = 0, 2 and w = 0. 

We now obtain for x(<7) the expression calculated to second order in 




(15) 


X(?)^(^)(4v)‘^|r«>(i)./[r«>(i)-i^[x>'^r«>(x)].j|^[x*^“r<§>(x)].^(|r5[x^'*r‘g>(x)]^^|. 

In view of the solutions obtained above, we finally express in the form 

xW. (=^) r (i^l (-|[x-'‘fMl.| [-|(x‘'«f wvj-fci) 




(16) 


This is the result of an exact solution to the irreducible equation, Eq. (6). In Sec. in, we will merely give 
the variational result calculated to the same order, and compare it with Eq. (16). 

ni. VARIATIONAL RESULT AND COMPARISON WITH OTHER RESULTS 

In Ref, 3, the vertex equation, Eq. (6), is solved by a variational method, and x((?) is computed for gen- 
eral wave vector q, frequency co. In this, we set co=0 and compute the static expression up to order. 
Without giving the details of computation, we quote the final result, here, in the present notation, 




(17) 


It may be of interest to point out that the variational 
solution to the vertex equation has the structure 

T.^{x;y)^i(y)Ai{y)-J(y)]. (18) 

where 

■^(y) = J dx Jdx 3^ (x;y) (18a) 

and 

J(y) = (ary2Tr^) j" dx ^ dx ^ x\^^ dxi j dxi 

X T(x, y)^{xi, y) 73( I x - Xi 1 ) 


I 

Here 

2i = l/[l+ (2ar./7r)iK«i>(l;l)] . (20) 

In deriving Eq. (19), one malces use of many sim- 
plifying properties that TM enjoys, even though 
the actual computation is somewhat tedious. 

In Sec. rv, we will discuss in detail the various 
aspects of a very popular model for Y^(lx-Xil), 
namely the Yukawa potential, 

IV. COMPUTATION OF gSeV OF THE GRADIENT 
EXPANSION 


X [1 - {xjx)^^^ cos^^^y cose,] . (18b) 

Observe that the trial solution rTir(^j‘3^) ^ 

dependence at all. When we compare Eqs. (16) 
and (17), we will be surprised to find that the dif- 
ference between them has a simple form: 


xW.x..(,Ux’(2^)lr-(5f.) 2. 

X ^ . (19) 


Shain^° has shown that If one writes x(.<l) in the 
form (we use the notations as in Sham’s paper) 

(21) 

then the coefficient of the gradient expansion 

in the theory of the inhomogeneous electron sys- 
tem is 

= [x‘®V(x^°’)* - Xa* V(Xo®’)®] , (22) 

where the subscripts 0 indicate the same quantities 
for the noninteracting system with the same density. 


GRADIENT EXPANSION IN THE DENSITY FUNCTIONAL.. 


If E„[«]i8 the exchange and correlation energy 
functional, Klelnman defined" 

^»[n| = J<„(«)ndV+ •• • , 

where €„(«) is the exchange and correlation energy 
per electron in the uniform gas. This serves to 
define The latter are easily calculated be- 
cause in the expressions for H we set 

the terms cont^nlng yJ" zero, we obtain 
respectively. 

We first make the observation that if we employ 
the first iterative solution to the T{fc;y) equation as 
was done by Sham, we obtain 

r(}ir;y)a 1 + j" lUi Y,(|x- Xi | ) 

Using this in Eq. (6) and the notations (18a, fSb), 
we obtain 

Xshain(<7)^ (^^jre^A7r^)[/(y) + J(y )] . (24) 

Even though the variational solution is of the form 
(18), in contrast to the iterative solution (23), we 
observe that XrtrM Is found to be 

_j 



Sham’s procedure is equivalent to e3q)an(liiig xiq) 
up to first order in or,/ tt, and so expanding XrirW 
up to this order. We observe this to be the same 
expression as Shames, 

On the other hand, if we use XrttU) to compute it 
to order 17®, we obtain 

and Sham obtains 

(27) 

This essentially means that the vertex renormali- 
zation r [3(1- is set equal to 1 + J^V 
I®’ in Sham's procedure. Since setting 11=0 and 
r = 1 imply that we are dealing with the noninteract- 
ing system, we obtain for the following expres- 
sions (in our notation): 

^g’(var) = 4irVwjS!j,e®/"’'r 

x{2(/®>//®>)[(l - r)/r] + J«>/7®>}, (28a) 
and 




4^’ 

(2ar.\ 


\ ^ ) 


2 

x-1 


Now from Eq. (17), we may deduce 

r=[l-(^)[Y®>(l;l)-iF<"(l;l)]]'\ 


and 


5^ 7 3 


(28b) 


(29) 


Thus we see that the variational result is ob- 
tained when Zi is set equal to zero, and Sham’s 
answer is obtained when Z^ is set equal to zero, 
and r equal to 1 + and working to first 

power in aTg/n (high-density limit) in the exact 
expressions. The correction to the variational 
answer is thus seen to be of order {otrj'n)^ includ- 
ing powers of otTg/Ti in Eq. (28b) without those ap- 


pearing in r and Z^, as they are of magnitude unity 
plus corrections of order ar^/TT provided that the 
expressions multiplying them are finite in this limit. 
In Sec, V, we will specialize the above results In 
terms of a Yukawa form for 7^(1 k - kj ) and dis- 
cuss the expressions in various regions of validity 
of such a model. We may state here that by ex- 
pressing our results in terms of as above, we 
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have made a contact with Ferml-liquid-theory 
parameters, which are being experimentally mea- 
sured these days. 

V. YUKAWA INTERACTION MODEL AND DISCUSSION 
OF THE RESULTS 


We will consider the effective interaction poten- 
tial between the electrons to be of the Yukawa 
form, which enables us to perform all the calcu- 
lations explicitly in terms of one parameter, the 
screening parameter. ^ Thus, 




and so, 


(30) 


X,) 

21+1 



(31) 


where Q; is the usual associated Legendre function. 
W& then obtain 




and 





Several features of this model may be of value and 
so we record them here: 

(i) In the large ^ limit (i. e. , short-range model), 

the vertex equation can be solved exactly and one 
has to observe that the coupling strength is now 
<xr,/T!f which is held constant. In this limit we 
have =1, but “ 0(1/?*), 

and so the variational result coincides completely 
with the exact answer, as is to be expected. 

(ii) In the limit of Coulomb interactions, i. e. 

0, we observe that goes to zero like 





( 1 .( 0 ) 


v-1/8 



diverges Uke (aV2ir)®[Ln(4/$®)f so that the cor 
rection to the variational answer diverges like 
(Q!ry27r) ln(4/|^). In the strictly Coulomb limits 
om has therefore a divergence for a fixed r, (or 


A similar result is also obtainod by Klein- 
man, But Klcinman takes to be the Thomas- 
Fermi expression, i.c., and claims 

that now there is no divergence, as Hence, 
as long as ( is treated as an independent parameter 
not depending on the electron density, one has a 
divergence, while if U dei>enda on this diver- 
gence disappears. Now ^ - 0 Implies in the Thomas- 
Fermi model for screening that and so, the 
only inference possible is that in the high-density 
limit, the corrections to the variational answer are 
again zero. Also, we may iwint out that Sham’s 
result is not expected to be correct for low densi- 
ties, where (or A>-*0)* 

(iii) In Ref. 13, we have developed a self-con- 
sistent scheme to determine A** Tho argument 
runs as follows. We know that in the limit of zero 
Qf gives 5^. Wo have computed ii(q) using' 

a Yukawa interaction, ami this now expression for 
it should bo consistent with the al>ove general re- 
sult. Thus, we obtain^* 

^ (4Tr/A>^) limx{^/) 
a-o 


or 

ta 


(t)]} 


This may be rewritten in the form 
or. 


(35) 


(i^S./ll*UJc-Aii.ln(l.f4A*JJ. (38) 


This equation can be solved graphically to deter- 
mine ns a function of or^/iT. We observe that 
0 only if ar^Tr- 0, i, e, only in the very higli- 
density limit. Thus the Tlwmas-Fermi reUUOfl- 
ship is reliable only for very high density of the 
electron gas. Also, for we observe the 

right-hand side of Kq. (36) becomes proportiODtl 
to u(?fc/4)- Thus in the low -density limit, the 
screening, though similar to the Thomas-Ferml 
expression, is different by a factor i. This is 
not a reliable answer since one expects for very 
low densities zero screening. 

It is worth noting that n(A>) ^ -(i^/ir*);f2\ cal- 
culated using our formalism, vanishes in the two 
extreme limits of itky - 0 (low density) and 
(high density) in contrast with the result of Klcin- 
man. Kleinman obtained a constant value for 
^{kp) for nkp - 00 (high density). To check if this 
would obtain in our formalism too, if we took f 
to be the Thomas-Ferml value 4ar^/ff, we made 
a similar asymptotic calcultitlon from our expres- 
sion for r2(Av), and we found this to vanish. For 
TT^jr-oo, we have tho extreme high-density limit 
where Sham’s procedure^*^ is certainly valid, as is 
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Lso evident from our ftnalysU. 

0 for and we are in agreement ^th 

Im* For 0 (low density), we agree with 
leinman's result that 0-^0 again, whereas Sham’s 
36uU is divergent but inappropriate in this limit. 
i view of the various comments made on the 
ukawa model, it sRjears a calculation of 
ith a more realistic momentum-dependent inter- 
3tion potential should be made. In Fig. 1 we have 
isplayed our results. It is found that the variation- 
L results are in excellent agreement with the exact 
asults for both high and low densities, except the 
jsition and height of the maximum are slightly 
Utferent, We have therefore displayed only the 
cact results for the self -consistent and Thomas - 
ermi schemes for f. The other contrastii^ fea- 
ires with the results of Sham and Kleinman are 
rident from this figure. A numerical table of 
3 Ttkjr may' be obtained on request, from either 
f the authors. 


N IN THE DENSITY FUNCTIONAL, 


Professor Kleinman has informed one of us 
(A, K.R. ) that he has now Verified that his D also 
vanishes for in agreement with Sham and 

our result. The general shape of his curve of D 
vs nkp now is the same as ours but seems to have 
some small differences in detaiL We thank Pro- 
fessor Kleinman for informing us of his new finding 
and for pointing out a sign error in the definition 
of O in terms of We must mention that we have 
also come across a paper by Geldart et aZ. [Solid 
State Commun. 16^ 243 (1976)], where a rigorous 
analysis of the also made, and their con- 

clusions in certain limits seem to be in general 
agreement with ours. These authors do not give 
a graph of 0 vs irkjr. 
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